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ABSTRACT 


An  analytical,  study  was  conducted  to  determine  if  available  unsteady 
normal  force  and  moment  aerodynamic  test  data  could  be  used  in  exjunction 
with  existing  helicopter  rotor  aeroelastic  and  variable  inflow  analyses  to 
provide  a  method  for  predicting  the  stall  flutter  response  of  a  helicopter 
rotor  blade.  For  this  purpose,  incompressible  unsteady  aerodynamic  data 
for  an  NACA  0012  airfoil  executing  pure  sinusoidal  pitching  motions  were 
employed.  To  apply  such  data  under  rotor  blade  operating  conditions  where 
multiharmonic  motions  and  velocity  variations  exist,  the  data  were  expressed 
as  functions  of  instantaneous  section  angle  of  attack,  angular  velocity, 
and  angular  acceleration.  In  addition,  scaling  procedures  were  developed 
in  an  attempt  to  account  for  the  effects  of  compressibility. 

Limited  application  of  the  resulting  analysis  to  define  the  aeroelastic 
characteristics  of  several  blade  designs  showed  that  significant  self- 
excited  torsional  oscillations  of  the  stall  flutter  type  could,  in  fact, 
be  predicted  for  certain  combinations  of  flight  conditions  and  blade 
designs .  Correlation  studies ,  to  evaluate  the  ability  of  the  analysis  to  predict 
control  loads,  were  performed  with  CH-53A  maneuvering  flight  test  data  and 
with  level  flight  test  data  from  the  NH-3A  (S-61F) .  Although  the  analyses 
produced  self -excited  stall  flutter  oscillations,  the  degree  of  correla¬ 
tion  indicates  the  need  for  further  development.  In  level  flight,  the 
oscillations  occurred  only  at  rotor  thrusts  in  excess  of  those  observed 
in  flight  tests.  Oscillations  were  predicted  during  certain  maneuvers 
with  good  agreement  with  the  persistence  of  stall  flutter  and  the  azimuthal 
onset  of  the  oscillation.  Although  essential  to  the  prediction  of  stall 
flutter,  unsteady  aerodynamics  are  shown  to  result  in  a  degradation  of 
blade  bending  stress  prediction,  which  appears  to  be  a  significant  factor 
affecting  control  load  correlation. 

It  is  recommended  that  studies  be  undertaken  to  investigate  the 
validity  of  a  basic  assumption  of  this  analysis;  namely,  that  unsteady  aero¬ 
dynamic  characteristics  of  an  airfoil  subjected  to  a  nonsinusoidal  angle 
of  attack  variation  into  stall  can  be  synthesized  from  unsteady  data 
obtained  for  sinusoidal  motions.  In  addition,  attempts  should  be  made  to 
improve  the  accuracy  of  the  analytical,  method  described  in  this  report 
through  possible  semiempirical  modifications  to  the  unsteady  data  and 
scaling  techniques  as  well  as  possible  refinements  to  rotor  trimming 
and  inflow  modeling  procedures. 
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INTRODUCTION 


The  subject  Investigation  was  conducted  Jointly  by  Sikorsky  Aircraft 
and  United  Aircraft  Research  Laboratories  under  U,  S.  Any  Aviation  Materiel 
Laborator5.es  Contract  DAAJ02-68-C-00U8 .  The  objective  was  to  deteradne 
whether  existing  unsteady  serodynaaic  test  data  and  rotor  aeroelmstic 
analyses  can  accurately  predict  stall- flutter-induced  control  loads. 

Tbe  iuftyttmnn  speed  of  helicopters  is  currently  United  by  the 
occurrence  of  stall  on  the  retreating  blades  of  the  rotor.  Such  stall 
manifests  itself  in  a  performance  degradation  and,  more  importantly,  in  an 
increase  in  the  vibratory  loads  generated  by  the  rotor.  In  particular,  the 
buildup  in  control  system  loads  as  steady-state  stall  angles  are  exceeded 
on  the  rotor  appears  to  be  critical  and  is  believed  to  be  directly  associ¬ 
ated  with  the  occurrence  of  stall  flutter  of  the  blades  (see,  for  example. 
References  1  through  U).  Stall  flutter  is  defined  as  a  self-excited, 
single-degree -of- freedom  torsional  oscillation  of  an  elastic  blade  which 
can  occur  when  the  blade  is  operated  at  angles  above  stall.  Such  oscilla¬ 
tions  impose  large  moments  cn  the  rotor  control  system  and  are  often  the 
determining  factor  in  limiting  the  maximum  speed  of  a  helicopter. 

It  is  now  will  established  that  classical  rotor  theories  are  incap¬ 
able  of  predicting  overall  rotor  performance  characteristics  in  stall 
(References  5  through  9)*  let  alone  the  details  of  rotor  loading  required 
for  the  accurate  definition  of  critical  vibratory  loads.  Such  theories 
generally  assuse  that  tvo-dimeneional,  steady-state  airfoil  data  are 
applicable  and,  further,  that  the  inflow  velocities  induced  at  the  rotor 
blades  by  the  vortex  system  (or  wake)  of  the  rotor  are  constant  with  respect 
to  both  time  and  space.  In  an  attempt  to  provide  a  more  suitable  analysis 
for  predicting  self-excited  control  system  loads,  available  two-dimensional 
unsteady  airfoil  data  have  been  integrated  into  an  existing  aeroelastic 
variable  inflow  analysis  for  rotors.  This  modified  analysis  will  hereafter 
be  referred  to  as  the  Modified  UAC  Rctor  Analysis.  The  unsteady  data 
originally  considered  for  this  purpose  were  those  presented  in  References  10 
and  11  for  airfoils  executing  pure  sinusoidal  pitching  motion  (i.e.,  motion 
having  but  one  harmonic  or  frequency  component).  However,  as  discussed  in 
this  report,  only  the  data  for  Reference  10  were  ultimately  used  to  obtain 
the  results  presented.  Since  rotor  blades  do  not  execute  pure  sinusoidal 
motion  as  they  penetrate  the  stall  regime,  a  unique  feature  of  this  investi¬ 
gation  was  the  conversion  of  the  available  sinusoidal  results  to  a  more 
general  fora  in  which  the  local  section  normal  force  and  moment  are  ex¬ 
pressed  in  terms  of  instantaneous  section  angle  of  attack,  angular  velocity, 
angular  acceleration,  and  Mach  number.  A  conversion  of  this  type,  or  its 
equivalent,  is  required  before  section  fcrces  and  moments  associated  with 
typical  rotor  blade  angle  of  attack  variations  can  be  computed  in  a 
rational  manner. 

The  portion  of  the  study  conducted  by  the  Research  Laboratories  had 
as  its  principal  objectives:  (l)  the  generalization  of  the  available 
sinusoidal  unsteady  data  to  the  form  described  above,  (2)  the  integration 
of  the  generalized  data  into  the  existing  rotor  analysis,  and  (3)  the 
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application  of  the  resulting  analysis  to  investigate  the  effect  of  certain 
blade  design  and  operating  conditions  on  blade  stall  flutter  character¬ 
istics.  In  addition,  since  there  was  some  evidence  that  sustained  stall 
flutter  could  be  encountered  in  certain  aircraft  maneuvers,  a  final  objec¬ 
tive  of  the  Research  Laboratories  efforts  involved  the  modification  of  the 
existing  rotor  analysis  to  include  the  effect  of  prescribed  fuselage  motions 
on  blade  aerodynamic  and  dynamic  forces. 

Sikorsky  Aircraft  conducted  a  correlation  study  vhich  had  as  its 
objective  the  determination  of  the  accuracy  of  the  analyses  developed. 

The  correlation  study  included  comparison  of  analysis  with  NH-3A  (S-6lF) 
level  flight  and  CH-53A  maneuvering  flight  test  data.  CH-53A  correlation 
used  a  constant  inflow  model,  while  both  constant  and  variable  inflow  are 
used  in  the  HH-3A  analysis.  In  an  attempt  to  account  for  rotor-wing 
interference,  a  semiempirical  modification  to  the  NH-3A  inflow  model  was 
included.  Further  studies  of  the  mechanism  of  stall  flutter  were  initiated 
because  of  instability  in  regions  of  the  rotor  disc  where  the  blades 
appear  to  be  free  from  stall. 
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METHOD  OF  ANALYSIS 


The  results  presented  in  this  report  vere  obtained  using  the  Modified 
UAC  Rotor  Analysis.  Hiis  section  of  the  report  describes  the  technical 
features  of  the  analysis,  indicates  the  nature  of  the  modifications  made 
as  part  of  the  present  investigation,  and  summarizes  the  assumptions  made. 


MODIFIED  UAC  ROTOR  ANALYSIS 


A  simplified  block  diagram  of  the  Modified  UAC  Rotor  Analysis  is  shown 
in  Figure  1.  As  indicated,  there  are  two  principal  parts  to  the  overall 
analysis.  The  first  is  the  Blade  Response  Program,  where  the  motions  cf 
the  rotor  blade  are  computed  for  a  given  combination  of  blade  control 
angles,  fuselage  motions,  and  distribution  of  induced  velocity  over  the 
rotor  disc.  The  second  part,  namely,  the  Circulation  Program,  is  used  to 
compute  the  velocities  (rotor  inflow)  induced  at  the  rotor  by  its  vortex 
system  (rotor  wake).  Because  of  the  complexity  of  the  rotor  aeroelastic 
problem,  it  is  necessary  to  perform  an  iteration,  using  the  programs  as 
shown,  to  insure  that  the  input  and  output  quantities  of  each  are  reasonably 
self-consistent.  The  principal  modifications  to  the  basic  analysis 
incorporated  as  part  of  this  study  (indicated  in  Figure  1  by  the  dashed 
arrows)  include  (l)  use  of  unsteady,  two-dimensional  normal  force  and 
pitching  moment  aerodynamic  coefficients  in  certain  blade  section  opera¬ 
ting  regimes  and  (2)  modification  of  the  blade  equations  of  motion  to 
include  aerodynamic  and  dynamic  forces  associated  with  prescribed  time 
histories  of  fuselage  motion  and  blade  control  angles .  A  more  detailed 
description  of  the  two  individual  programs  is  presented  below. 


Blade  Response  Program 

This  program  determines  the  fully  coupled  response  of  a  flexible 
rotating  blade,  given  the  distribution  of  induced  velocities  over  the  disc, 
the  blade  control  angles,  and  the  fuselage  motions.  The  blade  deflection 
is  expanded  in  terms  of  its  uncoupled  natural  vibratory  modes  (normal 
modes).  The  normal  mode  technique  is  widely  used  to  solve  aeroelastic 
problems  (see,  for  example,  Reference  12,  p.  125)  and  facilitates  the 
numerical  integration  of  the  blade  equations  of  motion  through  elimination 
of  dynamic  coupling  terms.  In  this  investigation,  the  blade  response  is 
assumed  go  be  composed  of  a  rigid-body  flapping  mode,  three  elastic  flatwise 
modes,  a  rigid-body  lagging  mode,  two  elastic  edgewise  modes,  and  two 
elastic  torsional  modes.  The  basic  differential  equations  of  motion 
governing  the  response  of  each  of  the  blade  modes  are  given  in  Reference  13 
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far  rotor  blades  having  chordvlse  mass  unbalance  and  operating  in  steady, 
level  flight.  The  dynamic  and  aerodynamic  terms  associated  with  accelerated 
f: * «ht  of  the  aircraft  are  given  in  Reference  14  for  mass  balanced  blades. 
For  this  investigation,  these  terms  were  added  directly  to  the  equations  of 
motion  of  Reference  13.  (This,  in  effect,  assumes  that  terms  involving 
cross  products  of  the  nondlmensional  fuselage  accelerations  and  the  chord- 
vise  distance  between  the  center  of  gravity  and  elastic  axes  are  of  second 
order  and  can  be  neglected.)  For  convenience,  the  specific  modifications 
to  the  Important  equations  of  Reference  13  are  given  in  Appendix  I. 

The  aerodynamic  model  used  in  the  analyses  of  References  13  and  lb  was 
based  on  the  use  of  steady-state,  two-dimensional  airfoil  data.  In 
accordance  vLth  quasi-steady  theory  (Reference  12,  p.  279),  (1)  the  local 

section  angle  of  attack  was  defined  by  the  velocity  components  at  its 
three-quarter  chord  point  and  (2)  the  theoretic -\1  damping  moment  in  pitch 
vas  included.  In  defining  the  local  angle  of  attack,  the  inflow  at  the 
rotor  induced  by  the  vortex  system  representing  the  rotor  wake  was  included 
using  an  analysis  similar  to  that  of  Reference  15  •  Before  this  basic  type 
of  analysis  could  he  applied  to  the  stall  flutter  problem,  it  vas  necessary 
to  Incorporate  a  better  representation  of  the  effects  of  the  shed  wake 
vorticJty  (vak  vorticlty  arising  from  timewise  variation  of  blade  bound 
worticity)  at  all  angles  of  attack,  particularly  those  above  steady-state 
stall.  ¥0  accomplish  this,  an  attempt  vas  made  to  replace,  where  possible, 
steady-state  airfoil  data  by  unsteady  data.  The  latter  were  derived  from 
available  data  for  &  two-dimensional  airfoil  executing  forced,  pure 
sinusoidal  pitching  motion  and  operating  at  a  constant  air  velocity  (here¬ 
after  referred  to  as  sinusoidal  data).  Under  such  conditions,  the  reduced 
frequency  parameter  k  sc«/2u  is  constant  with  time.  Conversely,  typical 
rotor  blade  section  operating  conditions  involve  variations  in  k  due  to 
changes  in  both  the  velocity  of  the  section  U  and  the  effective  frequency 
of  the  motion  <*1  that  occur  as  the  blade  rotates .  A  rigorous  method  for 
applying  sinusoidal  data,  particularly  above  stall,  has  not  been 
established.  In  this  investigation,,  it  has  been  hypothesized  that  the 
sinusoidal  data  could  be  generalized,  through  cross  plots,  to  functions  cf 
section  instantaneous  angle  of  attack,  angular  velocity,  angular 
acceleration,  and  Mach  number.  Given  the  local  values  of  these  parameters 
for  each  section,  the  unsteady  lift  and  moment  of  each  section  can  be 
computed. 

As  originally  envisioned,  unsteady  normal  force  and  moment  data  from 
References  10  and  11  for  the  NACA  0012  airfoil  were  to  be  generalized  to 
functions  of  the  parameters  noted  above,  with  the  Mach  number  range  covered 
extending  up  to  0.6.  However,  it  was  found  that  the  test  matrix  employed 
in  Reference  11  vas  not  sufficiently  complete  to  permit  conversion  of  the 
data  by  the  process  described  in  Appendix  II.  Therefore,  an  attempt  was 


4 


made  to  scale  the  Inexpressible  data  of  Reference  10  to  approximate 
compressibility  effects.  The  scaling  relations  that  were  employed 
basically  assuae  that  the  ratio  of  actual  angle  of  attack  to  steady-state 
stall  angle  Is  the  critical  persaeter  that  aust  be  duplicated.  This  Is 
discussed  in  acre  detail  in  a  later  section.  Provision  vas  aade  In  the 
program  for  Investigating  the  sensitivity  of  the  results  to  certain  of  the 
parameters  introduced  by  the  scaling  procedures.  In  particular  the  Mach 
number  range  over  which  the  incompressible  unsteady  data  are  scale**.  The 
results  of  the  sensitivity  studies  are  also  discussed  in  a  later  section. 
Conventional,  steady-state  lift  and  moment  data  for  the  HACA  0012  airfoil 
section  ve*e  used  for  those  combinations  of  angle  of  attack  and  Kach  mrnber 
far  which  either  unsteady  data  were  not  available  or  the  scaling  procedures 
were  suspect.  Thus,  for  example,  steady-state  data  were  used  far  reverse 
flow  operating  conditions  because  applicable  unsteady  data  were  not 
available.  Also,  because  of  the  effect  of  Kach  number  on  airfoil  stall 
characteristics,  use  of  scaled,  incompressible  unsteady  data  at  conditions 
involving  Mach  numbers  above  0.5  and  angles  of  attack  above  steady-state 
stall  was  not  considered  valid.  Finally,  steady-state  drag  data  were  used 
throughout  the  program.  Further  discussion  of  the  treatment  of  the  unsteady 
data  is  given  in  the  section  of  the  report  entitled  UNSTEADY  AIRFOIL  DATA. 

The  modified  modal  equations  are  nonlinear  and  are  solved  using  a 
numerical  integration  technique  to  permit  realistic  coupling  between  the 
airloads  and  the  blade  response.  The  procedure  thus  has  the  capability  of 
analyzing  transient  as  well  as  forced-response  motions.  The  calculation 
generally  proceeds  by  alternately  computing  the  modal  accelerations  (from 
the  generalized  forces  exciting  each  mode)  at  a  given  blade  azimuth  angle 
f  and  numerically  integrating  these  accelerations  over  a  small  azimuth 
increment  A  ^  to  determine  the  modal  displacements  and  velocities  at  the 
next  blade  azimuth  (  ^  +  A^  ).  To  avoid  possible  mathematical  convergence 
problems  with  such  a  procedure,  it  is  necessary  that  (1)  A  ^  be  small  and 
(2)  the  accelerations  of  a  given  mode  not  be  expressed  in  terms  of  the 
acceleration  of  apy  other  mode.  These  requirements  present  no  difficulty 
with  a  normal  mode  type  of  analysis  when  the  rotor  is  in  steady  flight  and 
when  conventional,  steady  aerodynamic  characteristics  are  use  a  to  define 
blade  airloads.  When  accelerated  flight  is  considered  in  conjunction  with 
the  particular  unsteady  ,herodynanic  formulation  of  section  normal  force  and 
moment  discussed  previously,  the  generalized  force  exciting  each  mode 
becomes  a  function  of  the  various  modal  accelerations  as  well  as  a  function 
of  the  fuselage  accelerations.  In  determining  the  results  presented,  it 
was  assumed  that  (1)  the  fuselage  accelerations  were  prescribed  in  advance 
(e.g.,  from  flight  test  results)  and,  initially,  (2)  the  blade  section 
angular  acceleration  parameter  B  used  in  part  to  define  the  unsteady  aero¬ 
dynamic  coefficients  could  be  determined  from  numerical  differentiation  of 
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the  previously  generated  local  angle  of  attack  tie*  history.  The  former 
assumption  is  believed  to  be  justified  because  of  the  relatively  low 
frequency  of  the  fuselage  compared  to  that  of  the  blades.  The  latter 
assumption  Implies  that  the  calculated  value  of  B  will  lag  the  true  value. 
Die  results  of  initial  calculations  indicated  that  significant  values  of  8 
were  generated  only  when  a  high-frequency,  stall  flutter  type  of  torsional 
oscillation  occurred.  Ibis  suggested  the  possibility  of  approximating  B 
by  the  expression 
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where  the  second  portion  of  this  expression  is  the  component  representation  - 
of  total  section  angular  acceleration.  This  approach,  which  permitted  the 
azimuthal  lag  in  B  to  be  reduced  to  the  order  A ^  (approximately  3  deg), 
was  used  in  obtaining  the  results  presented  herein.  (Further  details  on 
the  concept  of  the  angular  acceleration  parameter  will  be  presented  in  the 
next  section  of  this  report.) 

Pertinent  output  from  the  Blade  Response  Progran  for  this  study 
consisted  primarily  of  integrated  rotor  loads  (thrust,  propulsive  force, 
etc.),  time  histories  of  section  operating  conditions  ( a  , a  , B  ),  section 
aerodynamic  coefficients  (C|,Cmc/4  »C<j  ),  blade  elastic  response,  and 
pitching  moment  at  the  blade  root  (i.e.,  pitching  moment  imposed  on  the 
control  system). 


Circulation  Program 

The  function  of  this  program  is  to  compute  the  circulation  distribution 
over  the  rotor,  given  the  distributions  of  section  angle  of  attack  a  , 
angular  velocity  parameter  A  ,  and  angular  acceleration  parameter  B  over  the 
blade,  as  veil  as  a  specified  geometry  of  the  rotor  wake.  For  this  study, 
the  wake  was  approximated  by  the  classical  skewed  helicoid  defined  when  the 
vortex  elements  generated  <’j  the  blade  are  convected  downstream  relative  to 
the  blade  at  a  resultant  ve.ocity  equal  to  the  vector  sum  of  the  blade 
rotational  velocity,  rotor  forward  velocity,  and  average  momentum  velocity 
through  the  disc.  Once  the  circulation  distribution  is  known,  the  inflow 
velocities  induced  at  the  blade  by  the  bound  and  wake  vorticity  of  the 
rotor  can  be  computed.  An  iteration  between  the  Circulation  and  Blade 
Response  Programs  is  then  used  to  assure  compatibility  of  the  inflow 
velocities  and  the  blade  aerodynamic  and  dynamic  boundary  conditions  (see, 
for  example.  Figure  1). 
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The  technical  approach  used  in  the  Circulation  Program  is  basically 
s dollar  to  that  described  in  Reference  15  and  generally  represents  a 
rotaiy-wing  equivalent  of  the  classical  lifting-line  approach  used  success¬ 
fully  far  fixed  wings .  There  are,  however,  some  differences  between  the 
UAC  Circulation  Program  and  that  described  in  Reference  15  •  These  are 
discussed  below. 

The  first  of  these  differences  is  the  elimination  of  all  shed  vorticlty 
elements  (elements  arising  from  variations  of  blade  bound  vortlclty  with 
time)  in  the  wake  model,  a  modification  that  contributed  substantially  to 
reducing  computing  time  without  significantly  altering  computed  circulations 
and  associated  inflow  velocities.  Thus,  the  only  vortex  elements  retained 
in  the  vake  arc  trailing  elements;  i.e.,  those  arising  from  spanwlse 
variations  in  bound  circulation.  As  in  Reference  15,  the  strength  of  the 
trailing  elements  is  permitted  to  vary  from  point  to  point  in  the  wake  to 
reflect  the  variation  in  bound  circulation  that  occurs  as  the  blade  rotates. 
It  is  believed  that  a  more  accurate  representation  of  shed  vortlclty  effects 
is  obtained  by  the  previously  mentioned  use  of  unsteady  airfoil  data  in  the 
Blade  Response  Program.  The  use  of  two-dimensional,  unsteady  data  basically 
implies  that  the  primary  effects  of  the  shed  vorticity  in  a  helical  rotor 
wake  are  due  to  the  wake  region  near  the  blade  and  thus  can  be  approximated 
by  those  in  a  fixed-wing  type  of  wake,  as  indicated  in  Figure  2.  Miller, 
in  Reference  1 6,  shows  that  this  approximation  Is  reasonable  at  rotor 
advance  ratios  ft  of  interest  to  this  study.  This  type  of  approach  greatly 
facilita+es  the  inclusion  of  nonlinear,  unsteady  aerodynamic  effects  due 
to  stall. 

Another  significant  difference  between  the  UAC  Circulation  Program  and 
that  of  Reference  15  lies  in  the  treatment  of  unsteady  effects  on  lift 
curve  slope  o,  angle  for  zero  lift  Oql,  and  stall  angle  .  In 

Reference  15,  steady-state  values  were  used  for  these  quantities.  In  the 
present  study,  o,  <*olj  811(1  amo*,u  were  considered  to  be  functions  of  the 
section  angular  velocity  and  angular  acceleration  parameters  as  well  as 
Mach  number.  As  in  Reference  15,  the  local  bound  circulation  is  assumed  to 
be  proportional  to  section  angle  of  attack  (measured  from  the  zero  lift 
angle)  until  am0XiU  is  reached,  following  which  no  further  increase  in 
circulation  is  permitted.  To  determine  the  functional  relationships 
involved,  the  normal  force  data  of  Reference  10  (also  presented  in  Table  I) 
were  plotted  and  approximated  as  in  Figure  3.  The  results  ere  shown  in 
Figure  k  for  incompressible  flow,  and  the  effect  of  Mach  number  on  these 
results  was  approximated  by  using  the  following  equations: 
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These  equations  state  that  (1)  the  angle  for  zero  lift  Is  Independent  of 
Mach  nosber,  (2)  the  variation  of  unsteady  lift  curve  slope  in  potential 
flow  with  Mach  number  follows  the  Prandtl-Glauert  relation,  and  (3)  the 
unsteady  stall  angle  of  attack  at  any  Mach  nuabex  can  be  determined  from 
the  assumption  that  the  variation  in  maximum  unsteady  lift  coefficient  with 
Mach  number  follows  the  same  variation  as  the  maximum  steady-state  lift 
coefficient. 

Finally,  provision  has  been  made  in  the  UAC  Circulation  Program  to 
eliminate  unrealistically  large  induced  velocities  which  would  otherwise 
be  computed  if  the  lifting  line  representing  the  blade  section  was  very 
near  a  vortex  element.  This  was  formally  accomplished  by  setting  such 
velocities  to  zero  whenever  the  distance  between  any  blade  section  and 
vortex  element  was  less  than  a  prescribed  distance  assumed  to  be  equal  to 
the  vortex  core  radius.  This  distance  was  selected  as  0.02R,  a  value 
believed  to  be  representative. 

It  is  recognized  that  factors  such  as  self -induced  rotor  wake  distor¬ 
tions  (Reference  17)  and  finite  blade  chord  effects  have  been  neglected. 
Their  inclusion,  however,  was  beyond  the  scope  of  this  investigation. 


SIMMARY  OF  ASSUMPTIONS 


Since  major  portions  of  the  Modified  UAC  Rotor  Analysis  used  in  the 
performance  of  this  study  are  documented  elsewhere  in  the  literature 
(References  13  through  15),  the  principal  assumptions  discussed  herein  and 
in  those  references  are  summarized  below  for  the  convenience  of  the  reader. 
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1.  The  aircraft  is  in  steady,  level  flight  or  in  accelerated  flight 
vhere  the  fuselage  velocities  and  accelerations  are  prescribed. 

The  rotor  rotational  velocity  is  constant. 

2.  Two-dimensional  unsteady  normal  force  and  moment  data  obtained 
for  an  airfoil  executing  forced  sinusoidal  pitch  oscillations 
under  constant  velocity  conditions  can  be  applied  under  rotor 
blade  operating  conditions  involving  velocity  variations  and 
multiharmonic  angle  of  attack  variations .  It  is  assumed  that 
this  application  can  be  accomplished  by  generalizing  the  basic 
data  to  functions  of  section  angle  of  attack,  angular  velocity, 
and  angular  acceleration. 

3.  The  effects  of  compressibility  on  unsteady  airfoil  characteris¬ 
tics  can  be  approximated  by  scaling  the  generalized  incompressible 
unsteady  data  using  scaling  procedures  that  generally  predict 

the  effect  of  compressibility  on  steady-state  airfoil  character¬ 
istics  . 

4.  Quasi-steady,  two-dimensional  aerodynamic  theory  is  applicable 
for  those  blade  sections  operating  at  conditions  of  high  Mach 
number  and  high  angle  of  attack  or  in  reversed  flow. 

5-  Section  drag  coefficients  are  given  by  steady-state  values. 

6.  The  velocities  induced  at  the  rotcr  by  the  vortex  system 

representing  the  rotor  wake  can  be  computed  using  a  lifting-line 
type  of  analysis  in  which  wake  distortions  are  neglected. 

7-  The  blade  has  an  elastic  axis  so  that  blade  deflections  can  be 
considered  as  the  superposition  of  two  orthogonal  translations 
of  this  axis  and  a  rotation  about  it. 

8.  Principal  blade  flexibility  effects  can  be  accounted  for  by 
considering  only  three  flatwise,  two  edgewise,  and  two  torsional 
vibratory  modes . 

9.  Blade  flap  and  lag  hinges  are  coincident  for  articulated  rotors. 

10.  The  local  center  of  gravity  is  assumed  to  lie  on  the  major 
principal  axis  of  the  section. 

11.  The  blade  is  linearly  twisted  along  its  span. 
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12.  Hie  folloving  quantities  can  be  assumed  to  be  I'mall  in  comparison 

to  unity: 

a.  Flap  and  lead  angles  (in  radians)  and  their  derivatives 

b.  Ratios  of  elastic  deflections  to  rotor  radius  and  their 
derivatives 

c.  Ratios  of  chordwise  distances  (i.e.,  chord,  center-of-gravity 
offset,  etc.)  to  rotor  radius 

d.  Built-in  twist  (in  radians) 

e.  Ratio  of  flap-lag  hinge  radial  distance  from  center  of 
rotation  to  rotor  radius 

f.  Froude  number  (q/Q  R  ) 

g.  Ratios  of  blade  thickness  dimensions  to  chord 

h.  Ratios  of  fuselage  angular  velocities  to  rotor  angular 
velocity;  ratios  of  fuselage  angular  accelerations  to 
square  of  rotor  angular  velocity. 

13.  On  the  basis  of  Assumptions  11  and  12,  the  following  types  of 

terms  in  the  equations  noted  can  be  neglected  as  higher  order: 

a.  Flatwise  and  edgewise  bending  equations: 

(1)  Second-order  products  involving  elastic  coordinates, 
fuselage  angular  velocities,  fuselage  accelerations, 
and  distance  between  blade  elastic  axis  and  center- 
of-gravity  axis . 

(2)  Ktird-order  products  involving  elastic  coordinates, 
chordwise  distances,  flap  angle,  lead  angle,  built-in 
twist,  and  flap-lag  hinge  offset. 

b.  Torsional  equation: 

(l)  Third-order  products  involving  elastic  coordinates, 

fuselage  angular  velocities,  fuselage  accelerations,  and 
distance  between  blade  elastic  axis  and  center-of- 
gravity  axis . 
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(2)  Fourth-ord.-r  products  involving  elastic  coordinates, 
chordvise  distances,  flep  angle,  lead  angle,  built-in 
twist,  and  flap-lag  hinge  offset. 

c.  Flap  angle  and  lead  angle  equations  -  second-order  tents 
involving  products  of  elastic  coordinates,  fuselage  angular 
velocities,  fuselage  accelerations,  ehordvise  distances,  and 
built-in  twist. 

d.  Section  velocity  equations: 

(1)  Second-order  products  of  elastic  coordinate.'. 

(2)  Third-order  products  involving  the  elastic  coordinates, 
chordwise  distances,  flap  angle,  lead  angle,  built-in 
twist,  and  flap-lag  hinge  offset  as  factors. 

All  equations  -  the  spamrt.se  component  of  acceleration  due 
to  gravity. 


e. 


UNSTEADY  AIRFOIL  DATA 


» 


DESCRIPTION  OF  ORIGINAL  UAC  DATA 


Hie  oscillating  airfoil  data  used  in  the  present  analysis  were 
obtained  in  an  experiment  previously  conducted  at  the  United  Aircraft 
Research  Laboratories  (UARL)  for  the  Sikorsky  Aircraft  Division  (Reference 
10).  Portions  of  these  data  were  used  ir.  a  previous  analysis  of  incipient 
stall  flutter  sponsored  by  USAAVLABS,  and  a  complete  description  of  the 
test  program  will  be  found  in  Reference  2.  A  brief  summary  of  the  test 
program  is  Included  herein  to  provide  necessary  background  for  subsequent 
sections  of  this  report. 

The  test  data  were  obtained  from  differential  pressure  transducers 
mounted  on  a  2-ft-chord  NACA  0012  model.  Hie  tests  were  conducted  in  the 
two-dimensional  channel  of  the  UARL  8-ft  main  wind  tunnel.  Data  were 
obtained  for  most,  but  not  all,  combinations  of  the  following  parameter 
values: 


Amplitude  of  motion,  a  =  -4,  6,  8  deg 
Mach  number,  M  =  0.2,  0.3 ,  0.4 

Mean  angle  of  attack,  aM  =  0  to  33  deg  in  3-deg  increments 
Frequency,  f  =0.0,  0.5,  1,  2,  4,  8,  12,  16  cps 

Hie  actual  points  for  which  data  were  taken  can  be  found  in  Tables  I,  II, 
and  HI  of  Reference  2,  for  M  =  0.2,  0.3,  and  0.4,  respectively.  It  can  be 
seen  from  these  tables  that  most  of  the  test  points  were  taken  at  the  lower 
Mach  numbers  and  at  the  higher  amplitudes.  This  fact  influenced  certain 
choices  that  were  made  in  the  course  of  the  data  analysis  and  which  will  be 
described  in  subsequent  sections . 


GENERALIZATION  OF  UAC  UNSTEADY  SINUSOIDAL  DATA 


Conversion  of  Data  to  <*  ,  A  ,  B  Dependency 

Earlier  in  this  report,  it  was  pointed  out  that  it  was  necessary  to 
operate  on  data  obtained  from  sinusoidal  tests  to  provide  a  set  of 
coefficients  applicable  to  operating  conditions  for  a  helicopter  rotor 
blade,  where  several  harmonics  of  motion  are  present.  No  definitive 
analytical  or  experimental  study  has  been  conducted  regarding  the  best 
procedures  for  applying  sinusoidal  data  under  such  conditions.  In  this 
investigation,  it  has  been  hypothesized  that  the  sinusoidal  data  can  be 
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generalized  by  converting  them  to  functions  of  the  instantaneous  angle  of 
attack  a,  angular  velocity  param<  jr  A ,  and  angular  acceleration  parametkr 
B  for  a  given  Mach  number .  A  necessary  (bat  not  a  sufficient)  condition 
for  the  validity  of  this  approach  lies  in  the  ability  to  characterize  the 
essential  features  of  the  sinusoidal  data  in  terms  of  only  these  three 
parameters .  While  such  is  the  case  in  potential  flow,  as  shown  below,  the 
ability  to  do  so  in  nonpotential  flow  had  not  been  completely  proven.  Even 
assuming  that  the  essential  features  of  the  sinusoidal  data  can,  in  fact, 
be  adequately  expressed  in  terms  of  a,  A  ,  and  B ,  it  remains  to  be  shown  that 
such  a  characterization  will  predict  unsteady  aerodynamic  characteristics 
in  nonpotential  flow  for  multi-harmonic  motions.  This  can  be  demonstrated 
only  by  correlation  with  experimental  results.  The  theoretical  background 
and  the  broad  concepts  of  the  data  conversion  technique  used  in  the  present 
study  are  described  below.  A  detailed  discussion  of  the  process  used  is 
given  in  Appendix  II. 


Consider  first  the  case  of  a  linearized,  potential,  incompressible 
flow  past  an  airfoil  oscillating  sinusoidally  in  pitch.  The  normal  force 
coefficient  Cn  and  the  pitching  moment  coefficient  C  m  can  be  expressed  in 
terms  of  three  constant  coefficients  which  multiply  the  angular  displacement 
a,  the  angular  velocity  a,  and  the  angular  acceleration  a  (see  p.  272  of 
Reference  12).  (The  theoretical  moment  coefficient  given  below  is  resolved 
about  an  axis  at  o0,  measured  in  semichords  aft  of  the  midchord  of  the  blade 
in  accordance  with  accepted  practice.  The  experimentally  determined  moment 
coefficient,  from  Reference  10,  is  referred  to  later  in  the  text  as  cmc/4  •) 


c 


* 

n 


2  7r  C  (k)a*+ 


l  +  2  C  (k) 


(M 


c»  =  *  (  ’■  5)  C  tk)  If  (  i  -  o0)|(v  i-)  CU)  -  £-j  <?-  f- 

In  these  equations  the  starred  quantities  are  complex, 
tradition,  the  Theodors en  circulation  function 


ter)V<)s“  (?) 

However,  following 


C  ( k)  =  F(k )  +  i G ( k ) 


(6) 


will  remain  unstarred. 
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The  complex  angle  of  attack  a*  may  he.  expressed  as  the  sum  of  the  mean 
angle  of  attack  <*m  and  the  instantaneous  displacement  from  the  mean  ae'wt 
or 


a*  =  aM'+5e'“'  (7) 

Ibis  expression  implies  that  the  Cn  and  Cm  responses  each  consist  of  a 
steady  part  dependent  on  the  mean  angular  position,  and  an  unsteady,  time- 
dependent  part  associated  vith  the  instantaneous  displacement  from  the  mean. 
These  steady-state  force  and  moment  expressions  are  not  included  in  Eqs. 

(4)  and  (5)  because  the  original  expressions  in  Reference  12  were  written 
for  displacements  relative  to  the  mean  position.  Hence,  in  making  use  of 
these  equations,  one  must  be  careful  to  include  an  additional  term  in  each 
of  them  to  account  for  the  steady-state  response.  This  is  done  later  in 
Appendix  ELI. 

Equation  (7)  has  a  real  part  given  by 

a  =  aM  +  a  coSwt  (8) 


and  the  first  two  time  derivatives  of  the  real  part  can  be  written  as 


—  oj Q  Sinud  = +um 


a  =  -u>2a  C0Swt  =  -<i>2(a-aM)  (10) 

The  second  form  of  each  equation  is  obtained  by  eliminating  the  time 
parameter  through  the  use  of  elementary  trigonometric  identities.  At  this 
point,  it  is  convenient  to  define  the  dimensionless  angular  velocity 
parameter  A  and  the  dimensionless  angular  acceleration  parameter  B,  which 
together  with  Eq.  (8)  make  up  the  fundamental  three-parameter  set  of  the 
present  study. 
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a 


aM  +  ffcoswt 


(ID 


A  =  ~ -  =  ±k  ya*-(  a  -  au)t  (IP) 

0  S  (Tu)2“  =  •  k2{0  *aM>  (13) 


vhere  k  is  the  reduced  frequency  for  sinusoidal  motion,  defined  by  the 
expression 


k  - 


C  u> 

2U 


(l*) 


Conceptually,  Eqs .  (k)  and  (5)  can  be  expressed  in  real  form  as 
functions  of  the  angular  displacement,  angular  velocity,  and  angular- 
acceleration;  hence,  both  Cn  and  Cm  can  be  functionally  written  in  terms  of 
the  fundamental  three-parameter  set,  a  ,  A ,  B ,  as 

cn  =  Cnlr,  A.B)  (15) 


Cm  =  Cm(  a,  A,  B) 


(16) 


The  conversion  of  Eqs.  ( k )  and  (5)  to  this  functional  form  foi'  potential 
flow  is  a  relatively  straightforward  task  and  is  done  in  Appendix  III.  In 
the  case  of  nonpotential  flow  (i.e.,  for  incidence  angles  greater  than  the 
stalling  angle),  the  relationships  are  no  longer  linear,  and  the  simple 
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expressions  of  Appendix  in  are  so  longer  valid.  However,  it  is  possible 
to  use  the  results  of  che  experiment  described  in  Reference  10  to  tabulate 
both  C»  and  Cmc/,  as  functions  of  these  three  parameters .  It  is  the  purpose 
of  this  section  and  cf  Appendix  n  to  show  how  such  a  tabulation  can  be 
constructed,  and  it  is  a  further  purpose  of  this  report  to  show  how  such  a 
tabulation  can  be  used  to  predict  stall  flutter  instabilities  of  helicopter 
rotor  blades. 

As  described  in  Appendix  H,  the  first  task  in  converting  the  sinus¬ 
oidal  data  to  a  generalized  fern  was  to  select  an  interpolation  parameter. 
The  frequency  f  was  found  to  be  the  most  convenient  parameter  for  this 
purpose,  and  a  series  of  cross  plots  was  then  prepared  for  C«nc/4  versus  f 
for  constant  values  of  the  remaining  parameters .  Specific  values  of  the 
fundamental  three-par aae ter  set, a ,  A ,  B ,  were  then  chosen  for  the  coordi¬ 
nate  points  of  the  final  tabulation.  The  values  chosen  were: 

a  =  0,  1,  2,  3....,  26,  27  deg 
A  =  0,  ^0 .01,  -0.025,  ±0.04 
8=0,  ±0.001,  ±0.004,  ±0.007,  ±0-01 

Otace  these  values  had  been  selected,  the  procedure  described  in  Appendix  H 
warn  followed,  in  which  data  points  were  Interpolated  versus  f ,  plotted 
versus  A,  faired,  cross -plotted  versus  B,  faired  again,  cross -plotted  versus 
a ,  and  faired  once  again.  Selected  points  from  these  smoothed  curves  were 
then  checked  against  comparable  points  from  the  original  sinusoidal  loops, 
and  adjustments  were  made  in  the  faired  curves  to  reduce  the  error  between 
the  original  loop6  and  the  final  faired  data.  The  agreement  between  the 
two  is  discussed  in  the  following  section. 

Finally,  these  unsteady  moment  coefficient  data  were  tabulated  and  are 
presented  in  Table  II.  Farther  details  of  the  tabulation  will  be  found  in 
Appendix  II.  The  original  table  of  unsteady  normal  force  coefficient,  first 
presented  in  Reference  10,  has  also  been  included  in  the  present  report  as 
Thble  I. 

A  few  selected  plots  of  Cmc/4  versus  a  are  shown  in  Figures  5,  6,  and 
7  for  B=  -0.01,  0,  and  0.01,  respectively,  and  the  three  curves  in  each 
figure  are  for  A  =  -0.025,  0,  and  0.025.  Values  of  Cm<;/4  for  zero  angle 
of  attack  were  obtained  from  the  theoretical  formula  discussed  in  Appendix 
HI.  In  each  of  the  figures  presented  here,  a  baric  behavior  or  trend  of 
the  curves  can  be  seen  (which  is  exhibited  by  the  entire  set  of  Cmc/4,  a 
cross  plots),  but  this  trend  is  most  clearly  observed  in  Figure  5  for  8  = 
-0.01.  The  most  striking  thing  about  the  dynamic  moment  response  is  the 
delay  in  the  occurrence  of  stall  when  the  angular  velocity  is  positive  and 
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the  reversion  to  a  stalled  condition  when  the  angular  velocity  is  negative. 
t*1k  is  clearly  seen  in  Figure  5  In  the  curve  for  A  0.025*  vhich  regains 
level  (and  hence  unstalled)  beyond  o=  20  deg  and  does  not  exhibit  any 
sharp  decreases  until  a=  25  deg.  These  values  of  a  are  considerably 
greater  than  the  steady-state  stalling  angle  of  a~  13  deg  (see  the  curve 
for  A  =  B  =  0  in  Figure  6) .  The  curve  for  A  =  0  lies  above  the  curve  for 
A  =  0.025  when  a  is  small  and  crosses  be’ov  it  vhen  a  is  large.  Similarly, 
the  curve  for  A  =  -0.025  lies  above  both  the  other  curves  for  small  a  and 
crosses  below  both  curves  for  large  a  .  A  similar  behavior  is  observed  in 
Figure  6  for  B  -  0,  although  the  maximum  angle  attained  before  the  curve 
departs  from  the  horizontal  (which  may  be  regarded  as  the  dynamic  stall 
angle)  is  smaller  for  each  value  of  A  than  the  comparable  value  attained 
in  Figure  5  for  negative  B  .  Finally,  Figure  1,  for  positive  8  *  shews  the 
same  typical  behavior,  although  the  curves  beyond  dynamic  stall  are  more 
irregular  than  those  shown  in  the  previous  two  figures.  Furthermore,  the 
dynamic  stall  angle  is  again  smaller  for  positive  B  than  it  was  for  either 
zero  B  or  negative  B  . 

This  graphical  inversion  of  the  order  of  the  curves  is  responsible  for 
the  loop  crossovers  (discussed  in  Reference  2  and  in  more  detail  in  the 
following  section  of  this  report),  which  contribute  to  the  negative  aero¬ 
dynamic  damping  in  stall  flutter.  To  illustrate  this  behavior,  refer  to 
Figure  8,  vhich  is  identical  to  Figure  7  except  that  the  right  halves  of 
two  hypothetical  moment  loops  have  been  superimposed  on  the  curves,  one  for 
the  potential  flow  regime  at  small  a  and  one  for  the  stalled  flow  regime  at 
large  a .  The  right  half  of  each  loop  has  been  chosen  because  this  is  the 
portion  for  vhich  the  angular  acceleration  is  negative .  Actually,  a  true 
loop  would  have  variations  in  B,  ranging  from  zero  at  the  mean  angle  of 
attack  to  maximum  negative  8  at  the  maximum  angle  of  attack.  However,  the 
use  of  hypothetical  loops  in  a  constant  B  plot  vill  be  permissible  for 
illustrative  purposes. 

First  consider  the  loop  in  the  potential  flow  regime,  limited  in  this 
case  to  a<  8  deg,  vith  numbered  points  1,  2,  and  3.  At  point  1,  the 
acceleration  has  just  become  negative  during  the  positive  A  portion  of  the 
motion.  At  point  2,  the  B  parameter  has  reached  its  negative  nuarimnm  value 
and  A  is  zero.  At  point  3>  B  is  again  approaching  zero  and  A  is  negative. 
Note  that  the  portion  of  the  loop  shown  here  is  elliptical  and  traverses 
the  CmCM  ,  a  plane  in  the  counterclockwise  direction.  As  shown  in 
Reference  2,  this  is  a  stable  loop  which  produces  positive  aerodynanic 
damping . 
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Now  consider  the  loop  in  the  stalled  flow  regime  with  numbered  points 
k,  5,  and  6.  The  description  from  the  previous  paragraph  concerning  the 
relative  values  of  A  and  B  can  be  applied  here  with  points  1,  2,  and  3 
replaced  by  points  4,  5,  and  6,  respectively.  However,  this  time  the 
portion  of  the  loop  shown  here  crosses  itself,  and  a  portion  of  the  loop  is 
traversed  in  the  clockwise  direction.  As  shown  in  Reference  2,  this  is  a 
destabilizing  portion  of  the  loop.  Under  some  circumstances,  this  portion 
can  grow  at  the  expense  of  the  stabilizing  portion,  and  the  result  is  a 
self-excited,  single-degree-of- freedom  stall  flutter  in  the  pitching  mode. 
However,  this  phenomenon  is  usually  self-limiting  because  enough  positive 
danping  can  be  produced  by  a  stable  left  half  of  the  loop  to  balance  the 
negative  damping  of  the  right  half  of  the  loop. 


Verification  of  Conversion  Process 


The  usefulness  of  the  Cmc/4  tabulation  hinges  on  its  ability  to  pro¬ 
duce  the  correct  value  of  the  unsteady  moment  coefficient  for  a  variety  of 
motions  typical  of  helicopter  rotor  blade  displacements .  A  necessary  (but 
not  sufficient)  condition,  then,  is  that  the  tabulated  values  must  produce 
reasonable  facsimiles  of  the  original  sinusoidal  loops.  As  part  of  the 
cro6S-plotting  and  smoothing  procedures  described  in  the  previous  section, 
a  number  of  loops  were  reproduced  from  the  tabulated  data  and  were  compared 
with  the  appropriate  original  loops .  Satisfactory  agreement  between  the 
original  and  reconstructed  loops  was  achieved,  and  the  results  are  shown  in 
Figures  9  through  13  for  a  variety  of  conditions. 

In  mr?t  cases  the  agreement  is  moderately  good;  in  a  few  cases  it  is 
excellent,  and  in  a  few  cases  it  is  poor.  However,  it  must  be  remembered 
that  each  original  loop  represents  only  a  single  cycle  of  an  oscillation 
that  exhibits  a  random,  stall-induced  characteristic  superimposed  on  the 
basic  wave  form.  Since  any  individual  reconstructed  loop  was  produced  by 
an  averaging  process,  it  is  not  surprising  that  there  are  individual 
departures  from  the  original  data.  Indeed,  it  is  encouraging  to  note  that 
there  is  an  overall  agreement  in  the  trends  exhibited  by  these  figures. 

Tnis  overall  agreement  is  further  amplified  in  the  discussion  of 
integrated  effects  which  follows  below.  First,  however,  it  is  instructive 
to  consider  the  individual  comparisons.  Figures  9  and  10  contain  loops  for 
constant  reduced  frequencies  of  k  =  0.112  and  0.225  (  f  =  4  and  8  cps), 
respectively,  which  are  frequencies  representative  of  the  middle  of  the 
original  test  range  and  are  probably  most  representative  of  the  helicopter 
operating  range.  As  stated  above,  the  agreement  is  generally  good,  par¬ 
ticularly  at  the  lowest  mean  angles  of  attack,  but  even  when  tla  curve 
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values  do  not  natch,  the  locations  of  the  curve  crossovers  are  in  reason¬ 
ably  good  agreement.  These  crossover  points  are  of  vital  importance  to 
the  stall  flutter  prediction  because,  as  described  in  Reference  10,  the 
transition  from  a  counterclockwise  enclosure  of  the  loop  to  a  clockwise 
enclosure  represents  the  transition  frca  positive  aerodynamic  pitch  damping 
to  negative  aerodynamic  pitch  damping. 

Figures  11,  12,  and  13  contain  loops  for  mean  angles  of  attack  of  15, 
13,  and  21  deg,  respectively,  and  indicate  the  ability  of  the  tabulation  to 
predict  the  effect  of  reduced  frequencies .  This  is  the  angle  of  attack 
region  which  is  most  susceptible  to  stall  flutter,  as  discussed  in 
Reference  2.  In  these  figures,  the  agreement  between  the  reconstructed 
loops  and  the  original  loops  is  generally  good  for  reduced  frequencies  of 
0.112,  0.225,  and  0.338.  The  agreement  is  only  fair  at  the  lowest  and 
highest  reduced  frequencies,  4  =  O.O56  and  0.450,  but  the  directions  and/or 
crossovers  of  the  loops  are  correct.  It  should  be  noted  that  in  Figure  13, 
for  aM  =  21  deg,  the  reconstructed  loops  for  the  highest  frequencies  exhibit 
flat  regions  in  the  vicinity  of  a  =  27  deg.  These  flat  spots  occur  because 
one  or  more  of  the  parameters  a  ,  A ,  or  B  have  exceeded  the  tabular  values 
given  in  Table  H,  and  the  computer  program  has  selected  "corner  values"  in 
these  regions.  The  parameter  values  at  which  these  corner  values  are 
selected  are  generally  outside  the  range  of  interest  in  helicopter  applica¬ 
tions,  and  these  discrepancies  are  not  believed  to  be  important  in  their 
effect  on  the  stall  flutter  prediction. 

As  a  final  check  on  the  overall  agreement  between  the  reconstructed 
loops  and  the  original  data,  the  aerodynamic  damping  parameter  Ha2  was 
calculated  for  the  reconstructed  loops.  Use  was  made  of  Eq.  (4o)  of 
Reference  2,  and  the  integrated  damping  for  the  reconstructed  loops  was 
compared  with  the  results  published  In  Figures  14  and  15  of  Reference  2  for 
a=8  deg.  This  comparison  is  prese  .ted  in  Figure  14.  The  results  shown 
here  corroborate  the  statement  made  earlier  that  the  overall  behavior  of 
the  reconstructed  loops  is  in  good  agreement  with  that  of  the  original  data. 
In  every  panel  of  Figure  14,  except  the  panel  for  the  highest  mean  angle 
( aM  =  27  deg),  the  agreement  between  the  two  sets  of  results  is  good.  In 
fact,  a  comparison  of  Figure  14  with  Figures  14  and  15  of  Reference  2  shows 
that  the  disagreement  between  the  damping  for  the  original  and  the 
reconstructed  loops  for  the  selected  amplitude  is  generally  less  than  the 
experimental  scatter  of  the  data  for  all  amplitudes  shown  in  Reference  2. 
(The  faired  curves  of  Reference  2,  representing  all  amplitudes,  are 
included  in  Figure  14  for  the  convenience  of  the  reader.)  Furthermore,  in 
the  last  panel  of  Figure  14,  for  aM  =  27  deg,  where  a  substantial  disagree¬ 
ment  is  observed,  the  damping  for  the  reconstructed  loops  is  in  agreement 
for  low  values  of  reduced  frequency  and  is  conservatively  destabilizing  for 
large  values  of  reduced  frequency. 
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Additionally,  some  comparisons  of  reconstructed  normal  force  loops  and 
the  original  loops  from  Reference  10  are  shown  in  Figures  15  and  16.  These 
figures  show,  respectively,  the  effect  of  varying  mean  angle  of  attack  at 
a  typical  value  of  reduced  frequency  and  the  effect  of  varying  reduced 
frequency  at  a  typical  mean  angle  of  attack.  As  with  the  comparison  of  the 
moment  loops,  the  agreement  is  generally  good,  although  individual  loops 
have  local  discrepancies .  One  significant  aspect  of  the  unsteady  normal 
force  at  angles  above  the  steady-state  stall  angle  is  the  tendency  to  over¬ 
shoot  tte  steady-state  maximum  normal  force  by  a  considerable  margin.  These 
figures  show  that  the  maximum  value  of  the  unsteady  normal  force  can  be  as 
large  as  2.2,  whereas  the  maximum  steady-state  value  at  lew  Mach  numbers  is 
approximately  1.2  to  1.3,  Tills  phenomenon  is  noted  later  in  the  text  in 
connection  with  observed  delays  in  rotor  lift  stall. 


SCALING  OF  UNSTEADY  DATA 

The  unsteady  aerodynamic  characteristics  of  an  HACA  0012  airfoil 
executing  forced,  pure  sinusoidal  pitching  motion  form  the  basis  of  this 
stall  flutter  study*  As  originally  planned,  experimental  data  describing 
the  two-dimensional  normal  force  and  pitching  moment  characteristics  over  a 
vide  subsonic  Mach  number  range  were  to  be  employed  to  provide  an  improved 
analytical  method  capable  of  predicting  stall  flutter.  The  unsteady  airfoil 
data  were  available  from  two  sources.  One  series  of  tests,  conducted  at 
UARL  a  number  of  years  ago,  is  described  in  Reference  2.  In  these  tests, 
most  of  the  data  were  obtained  at  Mach  numbers  less  than  0.4  and  at  reduced 
frequencies  less  than  0.45.  More  recent  tests,  conducted  by  Boeing  under 
USAAVLABS  sponsorship,  are  documented  in  Reference  11.  In  the  Boeing  tests, 
the  Mach  number  range  was  extended  to  0.6  and  the  reduced  frequency  range 
was  extended  to  0-72. 

As  described  earlier,  the  application  of  unsteady  data  to  rotor  blade 
operating  conditions  required  that  airfoil  characteristics  be  expressed  as 
functions  of  four  variables:  Mach  number  M,  instantaneous  angle  of  attack 
a,  angular  velocity  parameter  A,  and  angular  acceleration  parameter  B.  In 
processing  the  Boeing  data,  it  was  found  that  the  test  matrix  employed  by 
Boeing  was  not  sufficiently  complete  to  permit  conversion  of  the  data 
( through  cross  plots )  to  the  required  M  ,  a  , A  ,  B  form  needed  for  this 
study.  (It  should  be  stressed  that  the  Boeing  tests  and  the  present  study 
were  not  designed  to  be  compatible,  although  agreement  between  comparable 
UAC  and  Boeing  data  points  was  generally  found  to  be  reasonable,  especially 
when  differences  in  steady-state  stall  angles  were  considered.) 
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To  satisfy  the  need  for  higher  Mach  number  data  in  the  form  required 
by  the  present  method,  a  Mach  number  scaling  procedure  based  on  the  steady- 
state  stall  angle  vas  developed.  In  the  absence  of  unsteady  data  at  higher 
Mach  numbers,  the  steady-state  normal  force  and  pitching  moment  character¬ 
istics  of  the  NACA  0012  airfoil,  normalized  in  terms  of  stall  angle,  were 
used  to  evolve  separate  scaling  laws  for  normal  force  and  moment .  These 
steady-state  scaling  laws  were  then  assumed  to  apply  to  the  unsteady  blade 
characteristics  as  well.  A  discussion  of  this  application  and  examples  of 
scaled  loops  are  contained  in  the  following  sections. 


Normal  Force 


First  consider  the  steady-state  normal  force  characteristics,  as  shown 
in  Figure  17  for  0.2  <  M  <  0.7-  These  data  were  obtained  from  lift  and 
drag  data  presented  in  Reference  18.  For  the  purpose  of  this  discussion,  a 
normal  force  "stall"  angle  asn  is  defined  as  that  angle  at  which  the  normal 
force  characteristic  at  any  M  initially  departs  from  linearity.  Also,  a 
dimensionless  normal  force  stall  angle  parameter  <rn  j  defined  as  the  ratio 
of  the  angle  of  attack  a  to  the  normal  force  stall  an^  :  asn. 


=  a/a*n 


(17) 


When  Figure  17  is  examined,  three  significant  variations  with 
increasing  Mach  number  are  clearly  evident  in  the  potential  flow  regime 
(i.e.,  a  <  asn where,  as  noted  above,  aSn  is  the  angle  at  which  the  normal 
force  curve  initially  departs  from  linearity) .  These  variations  are  an 
increase  in  normal  force  curve  slope  in  the  potential  flow  regime  (a<  asn  )} 
a  reduction  in  the  value  of  C  n  at  the  stall  angle  asn  f  and  a  reduction  in 
the  normal  force  stall  angle  itself.  By  properly  manipulating  these  three 
factors,  the  higher  Mach  number  curves  can  be  derived,  in  the  potential 
flow  regions,  from  a  single  incompressible  normal  force  characteristics 
curve.  In  particular,  it  can  readily  be  shown  that  in  the  potential  flow 
region  the  steady-state  normal  force  at  any  Mach  number  can  be  related  to 
the  incompressible  normal  force  by  the  expression 


(CiO 


°n , M  »  o 


v/l  -  M2 


(C"  STAU..8«)m 
(Cn  STALL,  Ss)m*0 


(18) 
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Biis  same  expression  was  then  assumed  to  apply  in  the  unsteady  situation. 
Note  that  Eq.  (18)  can  also  he  written  in  the  following  alternate  fora: 


J\  -  M2 


(c„). 


VL 


(cJ, 


(CnST*u..*s)n  (CnsnkLL,s»)i«»o 


(19) 


'nils  expression  effectively  states  that  at  any  Mach  number,  the  ratio  of 
the  normal  force  coefficient  (steady  or  unsteady)  to  the  steady-state  normal 
force  coefficient  at  the  stall  single  asn  is  equal  to  the  ratio  at  M  =  0, 
provided  both  ratios  are  evaluated  at  the  same  value  of  o-n. 

As  a  first  limited  test  of  the  applic ability  of  this  scaling  technique, 
Eq.  (18)  was  used  in  conjunction  with  the  incompressible  Cn  data  of  Bible  I 
(  a  =  B  =  0)  to  compute  steady-state  lift  characteristics  at  several  Mach 
numbers.  These  were  then  compared  with  actual  data  given  in  Figure  17 . 

The  comparisons  are  shown  in  Figure  18,  where  the  values  of  asn  and 
(  cnST4LL  ss  )  used  in  scaling  are  also  noted.  Below  asn>  agreement  between 
measured  and  scaled  results  is  excellent,  as  might  be  expected.  Bie  agree¬ 
ment  above  asn  is  less  faithful,  but  acceptable. 

As  a  further  check  on  the  general  accuracy  of  the  scaling  techniques, 
several  normal  force  hysteresis  loops  were  generated  for  M  =  0.4  and  0.6. 
These  were  compared  with  experimental  results  from  Reference  11.  In 
obtaining  the  scaled  results,  the  lata  of  Table  I  were  used  in  conjunction 
with  Eq.  (l8)  and  values  of  asn  and  (  CnST*LLiS,  )  for  the  Boeing  airfoil 
(see  Table  III).  Figure  19  compares  loops  for  M  =  0.4,  aM~  9.8  deg,  and 
a  ~  4.8  deg.  Experimental  loops  are  shown  as  solid  curves,  while  the 
dashed  curves  represent  loops  constructed  and  scaled  from  the  UAC  incom¬ 
pressible  unsteady  normal  force  data  of  Table  I.  Figure  20  gives  similar 
results  for  M  =  0.6.  Agreement  between  experimental  and  scaled  loops  is 
seen  to  be  reasonable  for  M=  0.4.  However,  at  M  =  0.6,  the  comparisons  are 
less  favorable. 

In  applying  the  scaling  techniques  for  the  investigation  reported 
herein,  the  following  computation  procedure  was  required  to  make  the 
transition  from  the  real  blade  conditions  to  the  data  tabulation  and  then 
to  scale  the  data  back  to  the  blade  conditions. 

1.  Compute  the  following  two-dimensional  conditions  at  a  rotor 
blade  section:  a  ,  M  ,  A  ,  B 

2.  Compute  the  dimensionless  stall  angle  parameter  an  =  a/asn 
where  asn  corresponds  to  the  blade  section  Mach  number  M 
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3*  Compute  the  equivalent  incompressible  angle  of  attack 
aMtO  =  (  asn^M=0 

4.  Enter  the  unsteady  normal  force  tabulation  with  aM.0  ,  A, 

B,  and  extract  the  normal  force  coefficient  M.0  for 

use  in  Eq.  (13) 

5.  Compute  a  maximum  normal  force  coefficient  at  stall 

(  cnST4LL  s,  )M  that  corresponds  to  the  blade  section 
Mach  number  m 

6.  Apply  Eq.  (13)  to  scale  (Cn)a  M.0  to  the  blade  section 
Mach  number  M 

7.  Result:  (cn)<rn,M 

■Hie  normal  force  coefficient  is  then  trigonometrically  converted  to 
the  lift  coefficient  at  each  spanwise  station.  An  integration  over  the 
blade  span  is  performed  to  determine  rotor  blade  loading,  and  the  sequence 
is  repeated  at  each  blade  asimuth  angle. 

The  scaling  parameters  used  in  applying  this  scaling  technique  to  the 
unsteady  data  in  the  blade  response  calculations  presented  herein  are  given 
in  Table  IV.  Note  that  for  M  =  0.2,  the  values  of  asn  and  (  cn$T4LLls  )  are 
different  from  those  used  in  the  calculation  of  the  M  =  0.2  curve  of  Figure 
18.  This  results  from  the  differences  in  steady-state  normal  force  response 
indicated  between  the  M  =  0.2  data  from  Reference  18  and  the  steady-state 
results  extracted  from  Table  I  for  A  =  B  =  0. 


Pitching  Moment 


As  in  the  case  of  the  normal  force,  the  scaling  procedure  used  for  the 
pitching  moment  is  based  on  the  steady-state  stall  angle.  Consider  the 
unpublished  steady-state  pitching  moment  characteristics.  Figure  21,  corre¬ 
sponding  to  the  lift  and  drag  data  of  Reference  18.  A  pitching  moment 
"stall"  angle  asm  is  defined  to  be  an  angle  near  the  point  where  the  steady- 
state  moment  at  any  M  departs  from  zero.  Also,  a  dimensionless  pitching 
moment  stall  angle  parameter  <rm  is  defined  as  the  ratio  of  the  angle  of 
attack  a  to  the  pitching  moment  stall  angle  asm.  Unlike  the  normal  force, 
the  pitching  moment  stall  angle  asm  is  not  necessarily  related  to  the 
departure  of  the  pitching  moment  from  linearity.  Rather,  it  was  initially 
defined  as  that  angle  which,  when  used,  permitted  the  steady-state  pitching 
moments  for  the  higher  Mach  numbers  (Figure  21)  to  be  predicted  with  reason¬ 
able  accuracy  from  the  incompressible  pitching  moment  data  (assumed  to  be 
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the  M  =  0.2  data  in  Figure  21).  Table  IV  lists  the  values  of  asm  found  to 
be  acceptable  on  the  basis  of  this  criterion  (see  original  aSm  column  in 
Table  IV).  The  actual  scaling  procedure  is  one  of  operating  on  the 
incompressible  data  with  the  following  simple  expression,  which  primarily 
manipulates  the  location  of  the  pitching  moment  breakpoint: 


(c"*|) 


■»<> 

\/i  -M* 


(20) 


Figure  22  compares  the  values  of  steady-state  moment  coefficients  calculated 
using  results  from  Table  IV  and  Eq.  (20) .  No  comparison  is  shown  for  M  = 
0.2,  as  the  agreement  is  exact  for  this  Mach  number.  Up  to  M  =  0.4,  the 
agreement  between  experimental  and  scaled  curves  is  seen  to  be  very  good. 

As  might  be  expected  from  previously  presented  lift  results,  the  agreement 
is  less  favorable  at  M  =  0.5  and  0.6. 


Hie  scaling  procedures  were  also  applied  to  generate  several  unsteady 
pitching  moment  hysteresis  loops  from  the  unsteady  data  of  Table  H  for 
comparison  with  experimental  loops  from  Reference  11  for  M=  0.4  and  0.6. 
Results  are  shown  in  Figures  23  and  24.  In  obtaining  these  results,  the 
values  of  asm  were  determined  from  the  Boeing  steady-state  moment  data 
(Reference  11)  and  are  given  in  Table  III.  The  comparisons  shown  in 
Figures  23  and  24  are  generally  not  favorable,  although  in  some  instances 
reasonable  agreement  in  loop  shape  and  crossover  points  can  be  noted. 

As  with  normal  force,  in  applying  the  steady  moment  data  of  Table  II 
in  the  Blade  Response  Program,  the  computation  procedure  must  first 
transform  from  the  real  blade  conditions  to  the  M  =  o  data  tabulation  and 
then  scale  back  to  real  blade  conditions .  A  description  of  this  sequence 
follows : 

1.  Compute  the  following  two-dimensional  conditions  at  a  rotor 
blade  section:  a  ,  M  ,  A  ,  B 

2.  Compute  the  dimensionless  stall  angle  parameter  crm  =  a/a^, 
where  aSm  corresponds  to  the  blade  section  Mach  number  M 
(see  Table  IV) 

3.  Compute  the  equivalent  incompressible  angle  of  attack 

aM=0  =  °m(a$m)*.o 
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4.  Enter  the  pitching  moment  tabulation  (Table  II)  with  0^  ,A, 

8,  and  extract  an  incompressible  moment  coefficient  (cmt/4)a 
which  appears  in  Eq .  (20) 

5.  Apply  Eq.  (20)  to  scale  (Cm*,^  Ml0  to  the  blade  section 
Mach  number 

6.  Result:  (Cm^)^  u 

m » 

This  sequence  is  repeated  for  adjacent  blade  sections  and  for  each  blade 
azimuth  angle. 

Following  the  initiation  of  this  investigation,  additional  data  on  the 
effect  of  compressibility  on  aerodynamic  damping  in  pitch  (equivalently, 
cmc/4)  vere  published  in  Reference  19-  These  data  are  presented  as  the 
dashed  curves  in  Figure  25,  where  the  boundaries  encompassing  the  negative 
damping  region  are  shown  as  functions  of  section  angle  of  attack  and  Mach 
number.  The  solid  curves  in  the  left-hand  panel  of  Figure  25  represent  the 
equivalent  boundary  defined  by  the  scaling  procedure  using  Cmc/4  data  from 
Table  II  and  the  original  values  of  from  Table  IV.  In  an  attempt  to 
improve  the  agreement  between  scaled  results  and  Reference  19  results,  a 
revised  aSm  tabulation  was  generated  (see  Tfeble  IV)  and  used  to  produce 
the  comparison  shown  in  the  right-hand  panel  of  Figure  25.  These  revised 
results  became  available  in  time  for  use  by  the  Sikorsky  Division  in  their 
correlation  studies.  However,  all  parametric  and  sensitivity  study  results 
presented  in  this  report  are  based  on  the  original  asm  tabulation  given  in 
Table  IV. 
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EXPLORATORY  STALL  FLUTTER  STUDIES 


« 


Tills  section  of  the  report  describes  the  results  of  exploratory 
calculations  made  to  assess  the  effect  of  several  blade  design  and  flight 
condition  persaeters  on  the  susceptibility  of  rotor  blades  to  stall 
flutter.  Presents  are  (1)  a  description  of  the  blade  designs  and  flight 
conditions  selected  for  investigation,  (2)  results  of  a  study  conducted 
for  a  reference  blade  design  to  establish  the  sensitivity  of  the  predicted 
torsional  response  (and  associated  root  .  tching  moment)  to  certain 
parameters  in  the  general  aerodynamic  model  used,  and  (3)  results  showing 
how  the  stall  flutter  characteristics  of  the  reference  design  are  affected 
by  changes  in  certain  design  parameters. 

^WDITIOMS 

The  parameters  expected  to  be  pertinent  to  the  stall  flutter  problem 
of  rotor  blades  are  listed  below: 

1.  Rotor  lift,  CL/ cr 

2.  Rotor  advance  ratio,  fj. 

3.  Rotor  propulsive  force,  Cpp/cr 

4.  Blade  twist  rate, 

5.  Blade  airfoil  section 

6.  Blade  elastic  axis  location 

?•  Blade  torsional  stiffness  (or,  equivalently,  torsional  frequency 
ratio,  ) 

8.  Blade  center-of -gravity  axis  location  (eg). 

Items  1  through  4  principally  influence  the  extent  to  which  the  steady- 
state  stall  angles  of  attack  will  be  exceeded  over  the  rotor  disc  for  a 
given  flight  condition.  Items  5  and  6  affect  the  unsteady  aerodynamic 
hysteresis  characteristics  of  the  blade  as  it  oscillates  above  and  below 
the  steady-state  stall  angle,  while  Items  7  and  8  contribute  to  the 
torsional  dynamic  response  characteristics  of  the  blade. 


BLADE  DESIGNS  ABD  FLIGBF 
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The  parameters  considered  in  this  investigation  were  Items  I,  k,  7, 
and  8  (lift,  twist,  torsional  frequency,  and  section  eg  location).  Items 
5  and  6  were  not  considered  because  of  a  lack  of  appropriate  unsteady 
airfoil  data;  all  calculations  herein  are  based  on  data  for  an  HACA  0012 
airfoil  oscillating  sinusoidally  in  pitch  about  its  quarter-chord.  In 
addition,  rotor  advance  ratio  and  propulsive  force  were  not  considered  as 
variables.  Rather,  as  discos sed  below,  typical  values  of  these  parameters 
at  which  the  stall  flutter  phenomenon  might  be  expected  to  occur  were 
selected.  Table  V  indi'^tes  the  blade  designs  investigated.  The  Reference 
Blade  was  substantially  an  S-61F  blade  (Reference  20)  except  for  its 
torsional  natural  frequency  (frequency  of  blade  when  mounted  on  control 
system),  which  was  increased  from  6. 5?  to  8.2P  (through  control  system 
stiffness  increases)  to  provide  a  greater  range  of  Ug  .  The  principal 
characteristics  of  this  blade  are  given  in  Table  VI.  All  other  blades 
studied  were  identical  to  the  Reference  Blade  except  for  the  parameters 
noted  in  Table  V.  In  the  case  of  the  Low-Frequency  Blade,  the  reduction 
in  frequency  was  achieved  by  reducing  the  control  system  spring  rate  from 
1.2  x  106  ft-lb/rad  to  0.07  x  10^  ft- lb/rad.  Figure  26  compares  the  mode 
shapes  and  frequencies  of  the  first  two  torsional  modes  for  these  two 
blades , 


All  blades  were  considered  to  form  parts  of  five -b laded  rotors  which 
were  operated  at  an  advance  ratio  fj.  of  0 .4  in  the  pure  helicopter  mode 
(i.e.,  rotor  providing  a  propulsive,  as  opposed  to  b  drag,  force).  A 
propelling  rotor  condition  was  selected  as  being  more  critical  from  the 
stall  flutter  standpoint  since  the  depth  of  rotor  stall  at  a  given  CL/cr 
increases  with  increasing  propulsive  force  (see  Reference  18).  Rotor  shaft 
angle  was  fixed  at  a  typical  value  for  the  propulsive  flight  mode  (  as~ 

-8  deg,  CpF  / or  ~  0.009  for  Reference  Blade),  and  the  torsional  response 
characteristics  of  each  blade  design  were  computed  at  rotor  lift  values 
corresponding  to  nominal  incipient,  moderate,  and  deep  s +all  operating 
conditions,  (Stall  here  is  defined  as  steady-state  stall.)  Increases  in 
lift  and,  hence,  penetration  into  stall  were  accomplished  by  means  of 
collective  pitch  (  675  )  variations.  'Jhe  required  values  of  075  were 
estimated  from  the  charts  of  Reference  18  for  the  Reference  Blade  as 


Condition 


e75,  deg 


Incipient  Stall  12.2 
Moderate  Stall  14.2 
Deep  Stall  15.2 


These  values  of  collective  pitch  were  also  found  to  be  appropriate  for  the 
other  blade  designs  considered.  Cyclic  pitch  control  angles  (  A,  and  t3<  ) 
were  arbitrarily  estimated  from  Reference  18  using  the  requirement  of  * 


27 


% 


approximately  zero  flapping  with  respect  to  the  shaft  for  the  Reference 
Blade  for  incipient  stall  operation.  The  cyclic  pitch  values  were  A,s  - 
-2  deg  and  =  9  deg. 

snBmvra  STUDIES 


As  originally  envisioned,  the  analysis  developed  in  this  investigation 
was  to  use  experimental  unsteady  data  for  Mach  numbers  up  to  0.6.  However, 
as  discussed  previously,  this  proved  to  be  impossible,  and  it  was 
necessary  to  use  scaling  procedures  in  an  attempt  to  account  for  the  effects 
of  compressibility.  Because  of  possible  uncertainties  introduced  by  such 
an  approach,  a  limited  number  of  calculations  were  made  for  the  Reference 
Blade  to  establish  the  general  nature  of  its  computed  torsional  response 
and  the  sen-  j  uivity  of  this  repons e  to  arbitrary  changes  in  certain 
parameters  introduced  into  the  aerodynamic  model  by  the  scaling  procedures. 
Also  included  in  this  study  w as  an  evaluation  of  the  sensitivity  of  the 
results  to  the  number  of  iterations  performed  between  the  Blade  Response 
Program  and  the  Circulation  Program  (see  Figure  1).  In  addition  to 
Indicating  the  relative  significance  of  the  various  parameters,  these 
calculations  also  aided  in  establishing  the  final  aerodynamic  model  used  to 
calculate  the  torsional  response  and  pitching  moment  characteristics  of 
the  remaining  blade  designs . 


General  Response  Characteristics 

Figure  27  indicates  (for  the  Reference  Blade)  the  general  effect  of 
Increasing  blade  collective  pitch  (equivalently,  rotor  loading)  on  the  time 
history  ol  blade  root  pitching  moment.  This  moment,  which  is  nondimension- 
alized  as  noted,  is  icposed  on  the  rotor  control  system  and  thus  is  a 
direct  measure  of  push  rod  or  pitch  link  loads .  To  a  good  first  approxi¬ 
mation,  the  moment  also  is  proportional  to  the  elastic  torsional  response 
of  the  blade.  For  comparative  purposes,  results  obtained  using  both 
conventional  steady-state  aerodynamic  data  (of  Reference  18)  and  scaled 
unsteady  data  (from  Tables  I  and  II)  are  presented.  Corresponding  time 
histories  of  angle  of  attack  at  the  three-quarter  radius  are  also  given  in 
Figure  28. 

The  results  of  Figures  27  and  28  indicate  that  increasing  075  causes 
the  angle  of  attack  on  the  retreating  blade  (  2J0  deg)  to  increase  and 

eventually  to  exceed  the  nominal  steady-state  stall  angle  (  ~  12  deg).  At 
the  lowest  ,  the  oscillatory  pitch  moment  is  small,  and  the  use  of 
unsteady  aerodynamic  data  has  little  effect  because  of  the  relatively  low 
section  angles  of  attack.  Increasing  collective  pitch  to  lb. 2  deg  results 
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in  an  oscillation  at  the  first  torsional  natural  frequency  of  the  blade 
when  unsteady  aerodynamic  data  are  used.  This  is  due  to  a  self- excited 
torsional  oscillation  of  the  blade  which  appears  to  be  initiated  at  about 
a  ^  of  220  deg.  The  oscillation  grows  in  amplitude  until  the  blade 
reenters  the  region  of  higher  dynamic  pressures  and  lower  angles  of  attack 
in  the  first  quadrant,  where  both  stall  and  unsteady  effects  are  reduced 
and  the  motion  is  quenched. 

The  torsional  oscillation  noted  is  characteristic  of  stall  flutter  and 
is  reflected  in  the  angle  of  attack  time  history  for  9n  -  Ik. 2  deg  in 
Figure  29.  The  oscillation  also  appears  to  be  superimposed  on  e.  lower 
frequency  moment  variation  which  is  partially  characterized  by  a  positive 
moment  peak  on  the  advancing  blade.  This  peak  results  from  a  relatively 
rapid  unstalling  of  the  blade  and  a  resulting  dynamic  overshoot  as  the 
blade  seeks  a  new  equilibrium  position. 

As  the  blade  collective  pitch  angle  is  increased  further  to  15.2  deg, 
the  characteristics  of  the  blade  moments  described  above  are  accentuated. 
Additional  information  about  the  self-excited  oscillation  noted  at  Bn  = 
13.2  deg  is  given  in  Figure  29,  where  the  relative  magnitude  and  phasing  of 
the  angle  of  attack,  angular  velocity  parameter,  angular  acceleration 
parameter,  aerodynamic  moment  coefficient,  and  elastic  twist  at  the  75 
percent  radius  station  of  the  blade  are  compared. 

When  steady-state  aerodynamics  are  used,  the  angles  of  attack  on  the 
blade  for  a  given  collective  pitch  value  are  higher  while  the  rotor  lift  Is 
lower  (see  Figure  28).  This  is  due  directly  to  the  occurrence  of  classical 
steady-state  stall;  however,  despite  the  apparent  deep  penetration  into  the 
stall  regime,  no  stall  flutter  type  of  oscillation  occurs.  The  higher  lift 
observed  with  unsteady  aerodynamics  is  of  some  significance,  in  that  this 
observation  is  qualitatively  consistent  with  delays  in  rotor  lift  stall 
noted  in  several  experimental  investigations  (e.g..  References  5  and  21). 

The  general  effect  of  including  the  velocities  induced  by  the  rotor 
wake  (variable  inflow)  in  the  unsteady  aerodynamic  model  is  shown  in 
Figure  30-  In  this  figure,  blade  pitching  moments  predicted  using  variable 
inflow  are  compared  with  corresponding  constant  inflow  results  from  Figure 
27.  This  comparison  indicates  that  the  blade  moments  predicted  with  vari¬ 
able  inflow,  while  exhibiting  some  differences  in  shape,  have  characteris¬ 
tics  essentially  similar  to  those  predicted  with  constant  inflow.  Stall- 
flutter-  type  oscillations  again  occur  at  fa  -  lk.2  and  15.2  deg  and 
have  the  same  frequency  and  approximately  the  same  amplitude.  Rotor  lift 
is  substantially  unchanged  by  the  inclusion  of  variable  inflow. 
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Effect  of  Variations  In  Aerodynamic  Model 


The  results  using  unsteady  aerodynamic  data,  presented  in  Figures  27 
through  30  and  discussed  in  the  preceding  section,  were  based  on  the 
following  assumptions  in  the  aerodynamic  model: 

1.  The  limit  Mach  numbers  (M(,  M, )  to  which  the  incompressible 
unsteady  data  could  be  scaled  at  angles  of  attack  above  and 
below  the  steady-state  stall  angle  were  0.5  and  O.85, 
respectively.  (Figure  31  indicates  the  parameters  used  to 
define  the  operating  regions  where  scaled  unsteady  aero¬ 
dynamic  data  were  used;  values  of  asn  and  asm  were  taken 
from  Table  IV.) 

2.  The  Prandtl-Glauert  compressibility  factor  Vl -Mz  could  be 
applied  at  all  blade  section  operating  conditions  regardless 
of  whether  they  involved  potential  or  nonpotential  flow. 

3.  All  results  using  variable  inflow  could  be  generated  using 
one  iteration  between  the  Blade  Response  and  Circulation 
Programs  (see  Figure  l).  This  was  equivalent  to  two  passes 
through  the  Blade  Response  Program  and  one  pass  through  the 
Circulation  Program. 

A  limited  number  of  calculations  were  made  to  examine  the  general  sensitiv¬ 
ity  of  the  blade  root  pitching  moment  to  each  of  these  assumptions .  Results 
for  toe  Important  14.2  collective  pitch  condition  are  presented  in  Figure 
32  and  are  discussed  below. 

The  effect  of  decreasing  M(  from  0.5  to  0.4  and  M2  from  O.85  to  0.6 
is  shown  in  the  top  panel  of  Figure  32.  The  results  indicate  that  both  the 
initiation  and  the  growth  of  the  self-excited  oscillation  in  the  third  and 
fourth  quadrants  are  substantially  unaffected  by  the  values  selected  for 
those  parameters .  This  results  frcm  the  faci,  that  the  local  angle  of 
attack -Mach  number  combinations  over  the  important  outer  region  of  the 
blade  are  such  that  unsteady  aerodynamic  characteristics  are  uied  for  both 
sets  of  M,  and  M2.  It  is  seen,  however,  that  in  the  first  quadrant  the 
selection  of  the  upper  values  of  V,  and  M2  leads  to  a  somewhat  larger 
high-frequency  torsional  oscillation.  This  is  due  to  the  more  extensive 
use  in  this  region  of  unsteady  data  with  their  inherent  destabilizing 
effects  at  high  angles  of  attack.  The  higher  values  for  the  scaling  limits, 
M,  =  0.5  and  M2  =  O.85,  were  arbitrarily  selected  for  use  in  all  remaining 
calculations . 
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The  effect  of  increasing  the  number  of  iterations  between  the  Blade 
Response  Program  and  the  Circulation  Program  (cf.  Figure  l)  from  1  to  2  is 
shown  in  the  center  panel  of  Figure  32.  As  indicated,  little  additional 
information  is  obtained  on  either  the  character  or  the  amplitude  of  the 
blade  moment  by  the  second  iteration.  A  single  iteration  was  therefore 
selected  as  being  cost  effective. 

The  bottom  panel  in  Figure  32  shows  the  effect  of  arbitrarily  removing 
the  Prandtl-Glauert  scaling  factor  for  those  conditions  where  the  angular 
velocity  parameter  A  of  the  blade  section  was  negative.  This  eliminates 
the  major  nonpotential  flow  operating  regimes  where  this  scaling  factor  is 
invalid.  As  shown  in  Figure  32,  the  moment  results  are  insensitive  to  the 
treatment  of  the  Prandtl-Glauert  factor.  For  the  remaining  calculations, 
the  Prandtl-Glauert  factor  was,  for  convenience,  applied  for  all  section 
operating  conditions . 

PARAMETRIC  STUDIES 

As  indicated  previously,  the  objective  of  the  parametric  studies  was 
to  investigate  the  effect  of  variations  in  blade  twist,  fivst-aode 
torsional  frequency  (i.e.,  torsional  stiffness),  and  chordwise  eg  position 
on  blade  stall  flutter  characteristics .  The  basic  results  of  the  calcula¬ 
tions  based  on  a  constant  inflow  assumption  are  presented  in  Figures  33 
through  35.  In  these  figures,  the  time  histories  of  blade  root  pitching 
moment  for  the  various  blade  designs  are  compared  with  those  for  the 
Reference  BLade  at  the  same  values  of  blade  pitcxi.  Note  that  rotor  lift 
for  a  given  pitch  setting  did  not  always  remain  constant  as  the  blade 
configuration  was  changed.  Therefore,  results  of  Figures  33  through  35 
have  been  cross-plotted  in  Figure  36  as  functions  of  rotor  lift.  In 
Figure  36,  the  |  peak-to-pesk  amplitude  (^PTP)  of  both  the  total  blade 
pitching  moment  and  that  portion  of  the  moment  corresponding  to  the  self- 
excited  torsional  oscillation  is  presented.  The  latter  amplitude  was 
arbitrarily  defined  as  the  maximum  amplitude  of  the  self-excited  response 
on  the  retreating  blade  (180  <  $  <  3&0  deg).  Larger  amplitudes  of  the 
self-excited  oscillation  are  noted  for  the  highest  pitch  setting  in  the 
first  quadrant;  however,  these  occur  in  a  region  where  the  results  are 
sensitive  to  the  choice  of  Mach  number  scaling  limits  (  M,  and  M2)  and  are, 
therefore,  su'  ject  to  some  question.  Results  similar  to  those  of  Figures 
33  through  36,  but  including  variable  inflow  effects,  are  presented  in 
Figures  37  through  40.  Discussion  of  all  results  follows. 
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Constant  Inflow  Results 


The  effect  of  blade  twist  on  blade  root  pitching  moment  characteris¬ 
tics  at  constant  pitch  setting  is  shown  in  Figure  33,  where  the  resets  for 
the  Reference  Blade  (  =  -4  deg)  are  compared  with  those  for  the  High 

Twist  Blade  (  6,  =  -8  deg).  (Hie  blades  were  identical  in  all  other 
respects . )  Increases  in  blade  twist  were  expected  to  be  beneficial  from  a 
stall  flutter  standpoint,  inasmuch  as  tLe  angles  of  attack  of  the  high- 
velocity  outboard  sections  cf  the  blade  are  reduced  at  a  given  rotor  lift. 
Close  examination  of  the  results  of  Figure  33  indicates  that  increasing 
blade  twist  does,  in  fact,  result  in  smaller  moment  oscillations  at  the 
same  rotor  lift.  Further,  the  reduction  is  primarily  due  to  a  delay  in  the 
onset  of  the  self -excited  oscillation.  Hiese  trends  are  more  clearly 
shown  in  Figure  36,  where  the  amplitudes  of  the  total  moment  and  self- 
excited  moment  component  are  presented  as  functions  of  rotor  lift. 

The  effect  of  varying  the  frequency  of  the  blade  first  torsional  mode 
(  "Og  )  is  shown  in  Figure  3^-  Results  are  presented  for  the  Reference 
Blade,  which  had  an  of  8.2P,  and  for  the  Low-Frequency  Blade,  which  bad 
an  u>0(  of  only  Up.  A3  noted  in  cn  earlier  section  of  the  report,  the 
reduction  in  wg  was  accomplished  by  arbitrarily  reducing  the  stiffness 
of  the  control  system.  This,  of  course,  altered  the  first  torsional  mode's 
shape  us  well  as  its  frequency  (see  Figure  26). 

Examination  of  the  results  for  the  Low-Frequency  Blade  leads  to  the 
following  important  observations .  First,  a  marked  reduction  in  the  rotor 
lift  far  a  constant  pitch  setting  occurs  when  Tig  is  lowered.  This  is 
due  to  the  larger  nose down  elastic  twist  angles  which  are  produced  by 
steady  centrifugal  and  aerodynamic  stall  moments .  Hie  magnitude  of  these 
moments  is  essentially  unaffected  by  wg(  and,  to  provide  equilibrium,  the 
less-stiff  blade  must  deflect  a  greater  amount.  A  second  noteworthy  point 
is  the  apparent  absence  of  large-amplitude,  self -excited  moment  oscillations 
that  could  be  clearly  attributed  to  stall  flutter.  Although  there  are  some 
moment  oscillations  at  a  frequency  approaching  that  of  the  first  torsional 
mode,  there  is  no  clear  indication  of  a  growth  and  decay  of  a  self-excited 
oscillation.  Hi  is  characteristic  appears  to  be  reasonable  if  one  assumes, 
on  the  basis  of  the  results  for  the  Reference  Blade,  that  one  or  two 
cycles  of  the  self-excited  oscillation  must  occur  before  significant 
amplitudes  are  developed.  Under  such  conditions,  a  self-excited  oscillation 
at  a  frequency  of  4P  may  never  have  an  opportunity  to  grow  significantly 
within  the  limited  azimuthal  range  where  negative  damping  exists .  Con¬ 
versely,  whatever  oscillation  does  develop  may  never  be  completely  quenched 
on  the  advancing  side  where  positive  damping  exists.  Hie  plots  of  moment 
amplitudes  as  a  function  of  lift  shown  in  Figure  36  Indicate  clearly  that 
the  Low-Frequency  blade  exhibits  significantly  lower  oscillatory  moments 


32 


at  a  given  lift  coefficient  than  does  the  Reference  Blade.  It  should  be 
remembered,  of  course,  that  the  comparison  would  be  less  favorable  if  made 
cn  the  basis  of  the  elastic  torsional  deflections  of  the  blades. 

Hie  final  design  parameter  investigated  was  the  blade  eg  chordvise 
position.  The  results  for  an  Aft  eg  Blade  (eg  at  30  percent  C )  are 
presented  in  Figure  3?,  where  they  are  compared  with  the  results  for  the 
Reference  Blade  (eg  at  25  percent  C  ).  All  moment  time  histories  for  the 
Aft  eg  Blade  are  characterized  by  a  particularly  severe  oscillation  on  the 
advancing  blade,  coupled  with  a  tendency  for  oscillations  near  the  first 
torsia  al  frequency  (  '  3P)  to  persist  over  much  of  the  azimuth  range. 

Bie  former  characteristic  is  a  steady-state  aerodynamic  phenomenon,  being 
observed  even  when  steady-state  aerodynamics  were  used  throughout  the  cal¬ 
culations.  It  appears  to  be  the  result  of  destabilizing  dynamic  forces 
associated  with  higher  harmonic  blade  flapping  and  bending  motions .  The 
time  histories  of  blade  pitching  moment  and  the  flapping  mode  acceleration 
are  ohcrwn  in  Figure  lt-1  for  both  the  Reference  Blade  and  the  Aft  eg  Blade. 
Bie  basic  characteristics  of  the  flapping  acceleration  are  the  same  for 
both  blades.  However,  for  the  Reference  Blade,  the  coupling  with  the 
torsional  mode  is  very  small  because  of  the  coincident  eg  and  elastic  axes . 
For  the  Aft  eg  Blade,  dynamic  coupling  exists  and  is  destabilizing  since 
the  torsional  response  produced  is  such  as  to  amplify  the  original  dynamic 
excitation.  The  high-frequency  oscillations  occurring  in  the  azimuth  range 
of  lte)  to  120  deg  are  not  due  to  stall  but,  rather,  are  believed  to  be 
transient  motions  induced  by  the  nearly  impulsive  dynamic  excitation  that 
occurs  on  the  advar.?irg  blade.  It  is  apparent  frem  the  results  of  Figure 
35  and  from  the  cross  plots  of  Figure  3 6  that  a  eg  axis  aft  of  the  blade 
elastic  axis  is  very  detrimental,  although  not  necessarily  from  a  stall 
flutter  point  of  vie/. 


Variable  Inflow  Results 

Blade  pitch  moment  results  obtained  when  variable  inflow  effects  are 
included  are  presented  in  Figures  37  through  40.  In  general,  the  previously- 
discussed  observations  and  conclusions  regarding  the  occurrence  of  stall 
flutter  for  the  Reference  Blade  and  the  relative  advantages  and  disadvan¬ 
tages  of  the  Lew-Frequency  and  Aft  eg  Blades  are  unaltered  by  the 
inclusion  of  variable  inflow.  However,  a  comparison  of  the  results  for  the 
Reference  Blade  and  the  High  Twist  Blade  is  somewhat  inconclusive.  The 
delay  in  the  onset  of  self-excited  moment  oscillations  with  increased  twist 
predicted  using  constant  inflow  theory  no  longer  appears  to  exist  when 
variable  inflow  is  included,  at  least  for  CL/cr's  up  to  0.082.  This  is 
probably  caused  by  the  unloading  of  the  blade  tip  region  due  to  the  tip 
vortex,  factor  that  makes  additional  unloading  associated  with  increased 
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twist  not  as  imnortant.  At  CL/<r's  greater  than  0.082  (e.g.,  for  *75  = 
15-2  deg),  the  computed  results  with  variable  inflow  indicated  a  marked 
reduction  in  the  amplitude  cf  the  self -excited  oscillation  for  the  High 
Twist  Blade.  Because  of  the  number  of  variables  involved,  the  exact  cause 
of  this  reduction  in  response  is  difficult  to  pinpoint.  It  was  noticed, 
however,  that  for  this  condition  the  phasing  of  the  aerodynamic  moments 
exciting  the  blade  tended  to  vary  more  along  the  span  of  the  blade,  thereby 
diminishing  their  effectiveness . 


CORRELATIOH  STUDIES 


The  objective  of  the  correlation  studies  was  to  determine  whether 
the  technology  available  in  unsteady  aerodynamics  and  rotor  dynamics  could 
be  used  to  accurately  predict  stall- flutter-induced  control  loads.  To 
evaluate  correlation,  two  significant  examples  of  stall  flutter  were  chosen: 
the  S-olF  (HH-3A)  in  high-speed  unaccelerated  flight  and  the  CH-53A  during 
maneuvers  at  high  load  factor.  The  stall  flutter  characteristics  in 
the  two  flight  regimes  differ  significantly.  The  S-61F  exhibited  strong 
oscillations  at  the  blade  torsional  natural  frequency,  which  initiated  over 
the  nose  of  the  aircraft  and  persisted  ever  the  retreating  side  of  the  disc 
(Figures  42-47 ) .  Attempts  to  duplicate  these  oscillations  analytically  in¬ 
cluded  examination  of  the  effects  of  airload  variations  (Figures  48-55)* 

In  addition,  the  effects  on  control  loads  correlation  of  collective  pitch 
setting,  unsteady  aerodynamics  and  blade  stress  predictions  are  considered 
in  Figures  56-60. 

The  CH-53A  maneuver  conditions  exhibited  oscillations  in  the  region 
of  the  torsional  natural  frequency,  which  persisted  over  a  major  portion  of 
the  disc  during  maneuvers  while  rotor  lift  and  control  load  amplitudes 
varied  (Figures  61-65).  In  some  cases  these  oscillations  persisted  fully 
around  the  rotor  disc  for  several,  revolutions  (Figure  66).  The  greater 
azimuthal  persistence  in  these  instances  is  attributed  to  the  high  rotor 
loading,  which  results  in  penetration  of  the  stall  region  by  a  significant 
portion  of  the  disc. 


S-61F  FLIGHT  TEST  CORRELATIOH 

The  S-61F  flight  test  program  provided  an  excellent  case  for 
correlation  because  of  the  clear  examples  of  several  cycles  of  high-fre¬ 
quency  oscillations  beginning  over  the  nose  of  the  aircraft.  In  addition, 
a  great  deal  of  in-flight  data,  including  airloads  data,  is  available  from 
the  USAAVLABS/HASC-sponsored  rotor  loads  measurement  program  (Reference  22) 
In  the  referenced  program,  the  airframe  and  rotor  system  were  heavily  in¬ 
strumented  to  obtain  in-flight  airloads,  hub  forces,  control  loads,  blade 
responses,  and  airframe  responses. 

In  performing  the  correlation,  a  variation  of  flight  speed  and  rotor 
thrust  was  considered  since  the  amplitude  and  azimuthal  persistence  of  the 
instability  are  dependent  upon  these  parameters.  The  conditions  analyzed 
are  referred  to  by  the  case  numbers  as  described  in  Reference  22.  Six 
cases  were  analyzed:  four  lS5-kt  cases  at  CT/cr  's  varying  from  .016  to 
.058,  and  two  165-kt  cases  at  C? /<r  's  of  .021  and  .054.  A  summary  of  the 
principal  parameters  describing  these  cases  is  provided  in  Table  VII. 

Mathematical  Model 

The  rotor  blade  was  modeled  with  two  rigid-body  degrees  of  freedom 
(flap  and  lag),  five  flatwise  elastic  modes,  two  edgewise  elastic  modes, 
and  two  torsional  elastic  modes.  The  modal  response  showed  that  the  energy 
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of  the  two  highest  flat  vise  modes  was  negligible,  and  results  were  not 
sensitive  to  the  deletion  of  the  fourth  and  fifth  flatwise  modes. 

When  constant  inflow  was  used,  analytical  rotor  forces  and  head 
acamnts  did  net  agree  with  flight  test  values  at  the  measured  settings  of 
cyclic  pitch,  collective  pitch,  and  shaft  single.  The  analysis  contains 
provisions  that  permit  the  computer  to  vary  pitch  in  a  prescribed  iterative 
manner  to  converge  on  lift,  propulsive  force  and/or  head  moments.  Of  the 
parameters  that  could  be  varied,  cyclic  pitch  values  were  the  most 
accurately  known  from  tests;  therefore,  the  sha^t  angle,  which  is 
extremely  sensitive  to  rotor-body  or  rotor-ving  interference,  vas  varied 
to  obtain  the  analytical  value  of  the  collective  pitch  required.  It  was 
determined  that  other  inflow  models  did  not  significantly  disturb  the  trim, 
end  the  settings  determined  using  constant  inflow  vere  retained  through¬ 
out  the  analysis. 

Experimental  and  Analytical  Push  Rod  Moments 

As  can  be  seen  in  Figures  42  through  47,  the  flight  test  push-rod 
mow*, at  increases  in  amplitude  with  increasing  rotor  lift  and  speed  and 
exhlMts  oscillations  at  a  frequency  in  the  vicinity  of  the  blade's 
torsiosai  natural  frequency  (approximately  6.5  eye les /revolution ) .  The 
oscillations  initiate  over  the  nose  and  continue  over  the  retreating  side 
of  the  rotor  disc.  Also  shown  are  analytical  moment  time  histories 
based  cn  a  constant  inflow  assumption,  which  do  rot  show  stall  flutter 
oscillation.  An  attempt  was  made  to  induce  the  oscillations  by  increas¬ 
ing  rotor  lift  beyond  the  test  values,  thus  increasing  the  penetration 
of  the  stall  region.  The  results.  Figure  56,  show  that  the  oscillation 
can  be  induced  in  this  manner,  although  the  instability  is  limited  to  the 
region  beyond  270  deg.  This  characteristic  is  similar  to  that  obtained 
in  the  sensitivity  study  and  again  substantiates  the  requirement  for 
unsteady  aerodynamics. 

Upon  noting  a  lack  of  correlation  vith  constant  inflow,  all  cases 
were  rerun  with  variable  inflow.  (The  variable  inflow  program  has  pre¬ 
viously  been  referred  to  as  the  Circulation  Program.)  Results  were  ex¬ 
amined  and  no  significant  change  in  the  character  of  the  response  was 
exhibited.  A  typical  comparison  of  constant  and  variable  inflow  results 
is  shown  in  Figure  43. 

The  presence  of  wings  and  auxiliary  propulsion  implies  rotor- 
ving  interference  and  smell  values  of  shaft  angle.  Rotor-ving  interference 
corrections  were  not  attempted  in  the  initial  correlation  studies  dis¬ 
cussed  above.  At  the  small  shaft  angles  characteristic  of  compound 
operation,  the  rotor  vakc  passes  close  to  the  disc.  In  such  cases,  the 
variable  inflow  program's  accuracy,  and  thus  the  resulting  airload  distri¬ 
bution,  is  somewhat  degraded. 

In  an  attempt  to  improve  the  correlation  between  predicted  and 
measured  airloads  for  the  S-61F,  a  combined  analytical  and  empirical  method 
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was  used  to  modify  the  rotor  inflow.  The  original  S-61F  correlation  cases 
were  run  with  constant  inflow.  Four  of  these  cases  were  rerun  with  vari¬ 
able  inflow,  using  the  same  inflow  analysis  as  in  the  sensitivity  studies. 

The  airload  distributions  resulting  from  these  cases  were  analyzed,  and  a 
further  modification  to  inflow  was  made  i.i  an  effort  to  match  the  test 
airloads.  The  modification  was  based  on  the  effect  of  wing  lift  on  rotor 
inflow.  The  inflow  due  to  the  wing  causes  a  downflow  over  the  rear  half 
of  the  rotor  disc  and  a  smaller  upflow  over  the  front  half.  This  inflow 
variation  was  approximated  at  each  blade  radial  station  by  an  azimuthal 
inflow  distribution  consisting  of  a  constant  and  one-  and  two-per-rev 
cosine  terms.  Relative  amplitudes  were  1 : 8 : U  for  the  steady,  one-per-rev, 
and  two-per-rev  coefficients.  All  three  terms  have  maximum  values  over 
the  tail  of  the  aircraft.  The  sum,  which  is  the  total  increment  on  rotor 
inflow  at  azimuth  zero,  was  chosen  to  provide  a  sufficiently  large  change 
in  blade  angle  of  attack  to  make  the  analytical  airloads  agr?e  with  the 
measured  airloads  at  the  tail.  This  modification  in  inflow  was  added  to 
the  previously  calculated  variable  inflow  distribution.  For  Case  U3, 
this  approach  had  to  be  modified.  For  this  flight  condition,  the  rotor 
wake,  as  given  by  the  variable  inflow  analysis,  remained  very  dose  to 
the  rotor  plane.  This  calculated  inflow  gave  large  distortions  in  rotor 
loads.  More  accurate  airloads  were  found  by  using  the  modification  method 
vith  the  constant  inflow  case  as  the  base.  The  only  harmonics  of  inflow 
used  for  this  case  were  the  one-  and  two-per-rev  components  required  to 
improve  airloads  at  the  tail. 

As  shown  in  Figures  U8  through  51,  this  procedure  did  improve  the 
airload  correlation.  The  change  in  blade  flapping  resulting  from  the 
inflow  modification  washed  out  part  of  the  change  in  angle  of  attack  so 
rotor  loads  at  the  tail  did  not  correlate  exactly.  The  airloads  in  the 
region  of  azimuth  90  deg  (over  the  right  wing)  did  not  show  much  improve¬ 
ment  in  correlation.  The  reason  is  probably  the  local  complexity  in 
inflow  due  to  the  proximity  of  the  wakes  of  the  preceding  blades  for  a 
compound  aircraft.  The  variable  inflow  analysis  used  does  not  consider 
wake  distortions  which  result  from  wake  interaction.  Thus,  the  true 
location  of  the  blade  wake  near  azimuth  90  deg  is  not  given  in  the  analysis . 
This  is  why  Case  1*3  could  not  be  used  with  variable  inflow.  Research  is 
under  way  at  United  Aircraft  for  the  development  of  a  program  which  does 
yield  wake  distortions.  This  will  permit  a  more  accurate  description  of 
airloads . 

The  cases  run  with  modified  airloads  showed  a  slight  increase  in 
the  oscillatory  push  rod  loads  but  did  not  improve  correlation  signif¬ 
icantly.  Results  sure  shown  in  Figures  52  through  55. 

Additional  Investigation  of  the  Stall  Flutter  Mechanism 

In  reviewing  the  results  of  the  Sensitivity,  Parametric,  and  Cor¬ 
relation  Studies,  a  number  of  facts  indicate  the  need  for  a  further  in¬ 
vestigation  of  the  stall  flutter  mechanism.  The  stall-flutter-type 
oscillation  was  achieved  in  the  Sensitivity  Study,  substantiating  the 
theory  that  unsteady  aerodynamics  are  necessary  for  stall  flutter  prediction. 
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However,  the  azimuthal  persistence  of  the  oscillation  is  significantly 
les«  than  that  exhibited  by  test .  It  vas  determined  that  the  use  of 
unsteady  aerodynamic  data  tends  to  reduce  angles  of  attack  for  the  same 
lift,  thus  reducing  the  penetration  of  the  stall  region  and,  consequently, 
control  loads.  This  phenomenon  is  exhibited  in  the  Correlation  Study, 
where  an  increase  in  rotor  lift  to  levels  above  test  was  needed  in  order 
to  initiate  the  oscillation.  An  attempt  to  modify  the  analytically  de¬ 
veloped  airload  distribution,  using  measured  airloads,  resulted  in  no 
significant  improvement  in  the  correlation.  Because  of  the  self-excited 
oscillations  observed  in  the  parametric  studies,  it  vas  concluded  that  the 
general  analytical  approach  can  lead  to  a  valid  prediction  method. 

However,  the  lack  of  substantial  aerodynamic  excitation,  with  unsteady 
aerodynamics  only,  implies  the  possibility  of  additional  effects  that 
cause  the  lack  of  correlation.  The  following  discussion  reviews  the 
results  of  additional  studies  into  the  stall  flutter  mechanism. 


A  major  consideration  is  the  inability  of  the  analysis  to  initiate 
the  oscillation  prior  to  azimuth  270  deg.  The  character  of  the  push-rod 
load  oscillation  is  somewhat  different  from  that  measured  in  flight.  The 
flight  test  oscillations  begin  earlier  in  the  rotor  disc,  before  azimuth 
180  deg  rather  than  in  the  vicinity  of  azimuth  270  deg. 

Since  the  blade  is  not  in  stall  on  the  advancing  half  of  the  rotor 
disc,  some  other  mechanism  appears  to  contribute  to  the  push-rod  load 
oscillation.  Torsional  moments  at  the  blade  root  can  be  generated  by  the 
coupling  of  flatwise  and  edgewise  response.  For  example,  among  the  com¬ 
ponents  of  torsional  moments  reacted  at  the  push  rods  are  those  produced 
by  the  product  of  blade  lift  and  edgewise  blade  displacement  due  to  bending 
and  by  the  product  of  drag  load  and  flatwise  displacement.  These  products 
are  referred  to  as  load-deflection  terms.  A  direct  determination  of  the 
importance  of  such  terms  could  not  be  made  from  flight  test  data,  since 
blade  deflections  were  not  measured  directly.  However,  blade  bending 
moments  are  a  measure  of  both  loading  on  the  blade  and  the  deflection  of 
the  blade.  A  product  of  the  flatwise  and  edgewise  bending  moments  at  a 
point  on  the  blade  serves  as  an  indicator  of  the  significance  of  the  load- 
deflection  terms.  If  the  azimuthal  variation  of  this  moment  product  cor¬ 
relates  with  the  push-rod  load  azimuthal  variation,  the  load-deflection 
effects  may  be  significant.  Figure  57  shows  a  typical  correlation  of  these 
terms.  Relative  magnitudes  have  no  importance  in  this  figure  since  the 
units  of  the  two  curves  differ  and  are  therefore  not  indicated.  However, 
it  is  significant  that  the  maxima  and  minima  of  the  moment  product 
consistently  lead  (by  a  small  phase  angle)  the  maxima  and  minima  of  the 
push-rod  load  curve.  In  the  region  from  270  deg  to  360  deg  the  moment 
product  has  a  small  value,  but  the  push-rod  load  oscillation  increases  in 
amplitude.  It  is  precisely  in  this  region  that  stall  flutter  oscillations 
were  predicted  by  the  analyses.  This  pattern  was  found  consistently  for 
the  flight  conditions  examined. 

This  correlation  of  moment  product  and  push-rod  loads  suggests  that 
the  oscillations  of  the  S-61F  push-rod  loads  are  initiated  before  the  blade 
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reaches*  the  stall  region  by  a  combination  of  flatwise  and  edgewise  response 
of  the  blade,  and  that  the  3tall  flutter  mechanism  further  amplifies  the 
push-rod  loads  on  the  retreating  half  of  the  rotor  disc. 

As  part  of  the  evaluation  of  the  correlation.  Case  43  was  rerun 
using  steady-state  aerodynamics.  Blade  pitch  and  aircraft  attitude  were 
held  the  same  as  in  the  unsteady  aerodynamic  analysis.  A  comparison  of 
the  two  cases  provides  same  insight  into  the  reasons  for  lack  of  correla¬ 
tion  of  push-rod  loads  (Figures  58  through  60/ .  As  compared  with  steady 
aerodynamic  results,  the  inclusion  of  the  unsteady  aerodynamic  terms  results 
in  greater  rotor  lift,  but  generally  lower  flatwise  and  edgewise  vibratory 
blade  stresses.  The  smaller  load-deflection  teiss,  which  result  from  the 
lower  stresses,  do  not  contribute  as  much  to  the  blade  torsional  response 
as  the  larger  load-deflection  terms  found  in  the  steady-state  case. 

k 

This  difference  is  evident  in  Figure  58.  which  shows  the  push-rod 
moment  for  both  the  steady  and  unsteady  aerodynamic  analyses.  The  use  of 
steady  aerodynamics  does  produce  a  small  oscillating  push-rod  moment  on 
the  retreating  half  of  the  disc,  while  no  such  oscillation  is  evident  with 
the  unsteady  aerodynamic  analysis.  It  was  shown  in  the  Sensitivity  Studies 
that  unsteady  aerodynamics  in  stall  will  induce  flutter.  That  no  such 
stall  flutter  occurred  with  unsteady  aerodynamics  in  Case  43  as  attributable 
to  the  low  angle  of  attack  on  the  retreating  blade,  so  that  the  unstable 
region  vas  not  sufficiently  penetrated  to  cause  flutter.  The  previously 
mentioned  load-deflection  terms  were  unable  to  act  as  the  oscillating 
mechanism  here  because  of  their  small  size,  although  they  had  done  so  in 
the  steady-state  aerodynamics  case,  where  they  were  larger.  One  may  then 
wonder  why  the  presence  of  unsteady  aerodynamics  only  was  sufficient  for 
flutter  in  Figure  28,  regardless  of  the  effects  of  the  load-deflection 
terms.  The  answer  would  appear  to  lie  in  the  higher  torsional  frequency 
used  there,  which  seems  to  desensitize  the  rotor  torsional  response  to 
loaa-aef lection  terms.  Figures  65  and  66  show  the  flatwise  and  edgewise 
stress  correlations  at  65  percent  and  46  percent  of  the  rotor  radius,  re¬ 
spectively.  Prediction  of  these  stresses  was  better  with  the  steady-state 
aerodynamics . 

Good  correlation  of  flatwise  and  edgewise  blade  stresses  apparently 
is  essential  for  good  correlation  of  push-rod  loads.  However,  even  when 
blade  bending  stresses  are  fairly  well  predicted,  as  in  the  steady  rero- 
dynamic  case,  something  more  is  needed  for  push-rod  load  correlations. 

From  the  Parametric  Studies  and  from  the  S-61F  cases  in  which  collective 
pitch  was  incremented,  it  is  evident  that  the  use  of  unsteady  aerodynamics 
can  produce  large  stall-flutter-induced  push-rod  moments.  This  effect  needs 
to  be  combined  with  an  accurate  prediction  of  flatwise  and  edgewise  moments 
to  provide  the  mechanism  for  initiating  the  oscillation  before  the  blade 
reaches  the  stall  region. 

It  is  not  evident  at  present  that  changes  need  to  be  made.  It  is 
possible  that  the  present  unsteady  aerodynamic  analysis  can  lead  to  correct 
results  if  rotor  trim  and  rotor  inflow  variations  are  considered  further. 

For  example,  the  poor  correlation  of  blade  stress  in  certain  azimuthal 
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regions  suggests  that  local  effects  such  cs  wake  distarticas  t  radial  flow, 
or  sirfimme-raScr  interference  can  be  most  important.  Possibly,  a  ciurs** 
is  needed  in  the  way  in  which  the  unsteady  aerodynamic  teres  are  used. 

This  nay  include  the  in  corporation  of  unsteady  drag  effects  or  a  modifica¬ 
tion  in  the  present  method  for  determining  blade  lift  and  pitching  moments. 

Tbc  results  of  the  S-6lf  correlation  show  that  the  analytical 
approach  being  used  is  generally  valid.  Continued  investigation  of  the 
area  of  rotor  trim  and  the  characteristics  «f  the  unsteady  aerodynamics 
of  ceciUaticg  airfoils  should  lead  to  an  accurate  analysis  for  the  predic¬ 
tion  of  push-rod  loads.  The  continuing  effort  should  include  tests  of  a 
two-dimensional  airfoil  executing  forced  pure  sinusoidal  motion,  forced 
miltihamonic  motion,  and  self-excited  stall  flutter.  These  tests  would 
he  designed  so  as  to  examine  the  validity  of  a  basic  assumption  made  in 
this  investigation;  namely,  that  unsteady  aerodynamic  data  for  pure  sinus¬ 
oidal  pitch  motions  can,  in  fact,  be  used  to  predict  -he  unsteady  char¬ 
acteristics  of  an  airfoil  executing  multibaroonic  pitch  motions  into  stall. 
In  addition,  the  tests  would  provide  detailed  stall  flutter  information  oc 
a  simple  model,  thereby  permitting  a  more  direct  evaluation  of  the  stall 
flutter  prediction  analysis. 

Additional  sensitivity  studies  should  be  undertaken  with  the  analyti¬ 
cal  method  developed  in  the  present  program.  The  objective  of  these  studies 
should  be  to  determine  if  the  unsteady  moment  data  of  Table  11  can  be 
modified  so  as  to  reduce  the  discrepancies  noted  without  destroying  the 
basic  ability  of  the  tabulation  to  reproduce  the  general  features  of  the 
original  sinusoidal  data  from  which  it  was  derived.  The  possibility  of 
developing  semiempirical  modifications  to  the  compressibility  ccaling  tech¬ 
niques  used  herein  should  also  be  examined.  Satisfactory  accomplishment  of 
this  objective  would  provide  a  stall  flutter  analysis  that  can  be  used  as 
an  interim  design  tool  while  more  fundamental,  long-range  investigations  of 
unsteady  aerodynamic  methods  are  pursued. 

MAHEUVERIH6  FLIGHT  TEST  CORRELATIOH  -  CH-53A 

The  CH-53A  flight  test  program  for  structural  substantiation  produced 
maneuvering  flight  data  under  conditions  that  had  produced  not  only  stall 
flutter  of  the  previously  observed  form,  but  also  cases  in  which  the  stall 
flutter  oscillations  persisted  over  much  larger  portions  of  the  azimuth. 

Two  cases  of  CH-53A  maneuvers  were  selected  for  correlation.  A 
short  transient  of  about  6  sec  duration,  involving  a  right  turn  at  a 
nominal  120-kt  speed  and  an  angle  of  bank  of  approximately  60  deg,  vas  the 
first  case  examined  because  it  could  be  considered  to  be  a  quasi-steady 
maneuver  for  most  of  its  length;  that  is,  although  it  was  accelerated 
flight,  the  accelerations  and  angular  accelerations  were  close  to  constant, 
and  there  was  little  variation  in  the  cyclic  and  collective  pitch  values. 
This  made  it  an  ideal  choice  for  initial  correlation  of  the  accelerated 
flight  provisions  of  the  analysis,  as  well  as  being  a  test  point  in  which 
stall  flutter  was  observed.  The  last  2  sec  of  this  flight  test  point 
involved  a  small  (approximately  3  to  5  deg)  increase  in  roll  angle, 
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rbicii  vsi  acccoyacied  by  eertre  vibrations  on  all  channels  of  data  ted 
extremely  high  control  ?oad  smplitudea .  However,  the  data  contained  too 
much  noise  or  vibrations  to  be  easily  used.  It  was  therefore  decided  to 
attespt  correlation  up  to  the  initiation  of  the  noise. 

The  second  wnaeuvet'  exa&ined  vas  a  longer  transient  of  sxproKite&ely 
23  sec  iavslviftg  a  «ominal  150-kt  right  turn.  The  angle  Of  bank  for 
this  point,  varied  fro®  under  30  deg  at  the  start  to  approximately  b0  deg, 
followed  by  a  rollout  back  to  level  flight .  The  more  complex  for*  of  this 
maneuver  served  to  provide  an  opportunity  for  observing  the  effects  of  u 
changing  load  factor  on  the  observed  stall  flutter  and  oa  the  accelerated 
flight  provisions  of  the  analysis . 

Hatheaatical  Model 

-  i  ■  i —  ■  i  W'*  n  *  •  i  i - r — - - - 

Ae  tr.  the  Sensitivity  and  Parametric  Studies,  the  modes  used  were 
tvo  rigid-body  degrees  of  freedom  (flap  and  lag)  and  three  flatwise,  two 
edgewise ,  and  two  torsional  vibratory  nodes. 

The  CH-53i  flight  test  aircraft  carried  a  large  number  of  acceler¬ 
ometers,  velocity  pickups,  and  inclinometers  at  various  locations.  After 
suppression  of  the  local  vibration  components  of  the  signals  and  other 
corrections ,  these  data  nere  processed  to  provide  16  time  histories  of 
motion  components  in  an  axis  system  affixed  to  the  rotor  hub.  These  time 
histories  represented  6  degrees  cf  freedom  of  aircraft  motion,  their 
necessary  time  derivatives,  and  histories  of  control  angles. 

The  maneuver  option  of  the  program  was  used  for  the  time  histories 
Just  mentioned;  that  is,  fuselage  motions  were  prescribed  along  with  time 
histories  of  pilot  hand  motions  on  the  cyclic  and  collective  controls.  The 
rotor  transient  response  was  examined,  rather  than  allowing  it  to  die  away 
tc  steady  state  as  in  other  sections  of  this  report. 

As  in  the  S-61F  cases,  flight-test-reported  cyclic  pitch  was  accepted, 
and  no  attempt  to  trim  head  moments  vas  made.  While  the  head  moments  of 
the  S-61F  were  measured,  those  on  the  CH-53A  were  unknown. 

The  CH-53A  had  carried  no  instrumentation  to  measure  rotor  lift, 
load  factor  was  available  either  from  measured  aircraft  weight  and  motions 
or  from  a  highly  damped  vertical  accelerometer  that  acted  as  &  load  factor 
readout.  Although  the  motion  and  load  factor  data  were  in  general  agreement, 
small  differences,  especially  at  the  higher  load  factors,  caused  acme  un¬ 
certainty  in  rotor  lift.  With  an  assumed  correction  for  fuselage  vertical 
drag,  the  rotor  lift  was  calculated  in  a  band  thst  varied  from  ±  3000  to 
±  6000  lb  about  the  mean  value  used  for  correlation.  This  represented  a 
band  of  about  ±8  percent  in  rotor  lift. 

Flight-test-reported  collective  pitch  and  shaft  angle  did  not  produce 
an  analytic  lift  near  the  calculated  experimental  rotor  lift.  Therefore, 
collective  pitch  vas  adjusted  until  rotor  lift  near  the  initiation  of 
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the  maneuvers  equaled  the  test  lift  as  determined  froa  gross  weicjot,  load 
factor,  aad  estimated  fuselage  vertical  drag,  Occe  this  was  accomplished, 
the  collective  pitch  was  varied  with  tias  in  a  manner  paralleling  that 
of  flight  teat;  that  is,  tv»  value  of  collective  pitch  reported  by  flight 
tost  at  the  instant  of  time  in  which  lift  wan  matched  was  net  used,  but  the 
departures  from  this  reference  point  were  followed  exactly. 

Only  constant  inflow  was  used  in  the  maneuvering  cases,  as  the 
present  variable  inflow  program  is  not  applicable  ,vhen  the  vortieity  pat¬ 
tern  is  tlae  varying v  Although  the  inflow  was  held  constant  over  the  disc. 
It  was  time  varying.  Co  the  basis  of  the  time  history  of  required  rotor 
lift,  itself  generated  from  the  load  factor  and  an  assumption  concerning 
body  vertical  drag,  a  momentum  value  of  inflow  was  computed. 

Dlscussl  on  of  CK-53A  Correlation 

As  in  the  S-61F  study,  the  degree  of  correlation  was  less  than 
satisfactory.  However,  good  correlation  was  achieved  in  portions  of  the 
maneuvers,  providing  aide  confidence  in  the  analysis.  Figures  6l,  63,  and 
64  show  the  variation  in  ^  peak-to-peak  push-rod  moment  and  rotor  lift  with 
time  in  the  two  transients.  In  the  short  (120-kt)  transient.  Figure  6l, 
the  analytical  rotor  lift  was  adjusted  in  an  attempt  to  match  test  values, 
which  are  shown  as  a  broad  band  because  of  the  uncertainty  of  the  fuselage 
vertical  drag  load.  Collective  pitch  and  rotor  shaft  angle  were  varied  in 
an  attempt  to  achieve  the  lift  calculated  to  have  been  present  in  flight. 
The  lifts  shown  are  predominantly  below  the  flight  test  band. 

Because  the  analysis  was  unable  to  generate  lifts  higher  than  those 
shevm,  it  i;*  probable  that  the  blade  is  heavily  stalled  over  much  of  the 
azimuth  and  has  penetrated  deeply  into,  or  possibly  through,  the  region  of 
negative  damping.  However,  the  fact  that  the  analytical  lift  is  lower  than 
the  flight  lift  would  imply  that  the  analytical  load-deflection  terms  would 
be  reduced.  Since  the  1/2  peak-to-peak  push-rod  moment  is  higher  for 
analysis  then  for  test  in  this  case  (Figure  6l),  it  can  he  inferred  that  the 
unsteady  aerodynamic  effects  predominate.  The  amplitude  of  this  moment  is 
largely  the  result  of  stall-flutter-type  oscillations  whose  character  is 
fairly  similar  to  that  of  flight  test,  as  can  be  seen  in  the  plot  of  a 
representative  revolution,  Figure  62. 

The  attempt  to  match  calculated  flight  lift  for  the  longer  (150-kt) 
transient  (Figure  63)  was  less  successful.  Although  the  correct  lift  was 
generated  in  revolutions  10  through  20  (the  area  selected  for  trim),  the 
remaining  portions  of  the  irtce  show  sizable  disagreements.  However,  the 
\  peak-to-peak  moment  trace  (Figure  64)  is  more  encouraging,  showing  trends 
similar  to  flight  tc-ot,  with  varying  degrees  of  amplitude  correlation. 

At  the  peak  valuer  (revolutions  40  to  50),  the  analytical  amplitude  is 
below  that  of  flight  test,  whereas  the  opposite  is  true  earlier  in  the 
maneuver. 

In  revolution  20  (Figure  65),  stall  flutter  oscillations  that  corre¬ 
late  well  with  flight  test  can  be  observed.  The  discrepancies,  which 
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account  for  the  higher  peak-to-peak  amplitude,  are  oc  the  advancing  side 
cf  the  disc.  They  sees  tc  reflect  slower  damping  oat  of  the  oscillation 
than  occurs  is  flight  test.  On  the  ether  hnnd,  revolution  «1  (Figure  66) 
shows  the  correct  amplitude  hut  withc  it  Jje  oscillation,  this  revolution 
occurs  iamediately  after  the  drop  in  lift  observed  in  Figure  63.  Such  a 
drop  would  indicate  both  a  lowered  angle  o?  attack  and  smaller  load-deflec¬ 
tion  terns.  Thus,  absence  of  oscillation  on  this  trace  is  not  surprising. 
The  cause  of  this  reduction  in  lift  has  not  been  established. 

In  general,  the  maneuvering  flight  correlations  emphasise  the  sane 
factors  as  the  S-61F  correlation.  To  correlate,  one  needs  both  unsteady 
aerodynamics  in  the  stalled  region  and  proper  load-deflection  terms 
to  initiate  oscillations  prior  to  stall.  Is  addition,  the  trend  of  the 
h  peak-to-peak  moments  in  Figure  61*  is  encouraging  and  tends  to  provide 
confidence  in  the  maneuver  modifications .  The  reasons  for  the  lift  dis¬ 
crepancies  (Figures  6l  and  63)  are  not  clear,  but  the  discrepancies  sees 
to  be  connected  with  the  constant  inflow  time  histories  used. 
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the  c one lu? iocs  listed  below  were  reached  in  the  course  of  this  study. 
Rote  that  conclusions  6  through  10  should  be  applied  only  in  the  determina¬ 
tion  of  possible  trends.  Further  validation  of  these  trends  mist  await 
analytical  improvements  leading  to  a  higher  degree  of  correlation  of  the 
analysis  with  flight  test  data. 

1.  Incompressible  unsteady  aerodynamic  normal  force;  and  moment 
characteristics  generated  by  an  HACA  0012  airfoil  executing  pure  sinusoidal 
pitch  motions  into  stall  can  be  expressed  with  reasonable  accuracy  as 
functions  of  instantaneous  angle  of  attack,  angular  velocity,  and  angular 
acceleration  (see,  for  example.  Figures  9  through  16). 

2.  Application  of  unsteady  data  in  the  form  described  in  (1)  above 
to  compute  rotor  blade  response  at  an  advance  ratio  of  0.4  leads  to  the 
prediction  of  significant  self -excited  torsional  oscillations  of  the  stall 
flutter  type.  These  oscillations,  which  occurred  above  a  rotor  lift 
coefficient/solidity  ratio  of  about  0.065,  directly  impose  large  oscillatory 
moments  on  the  control  system. 

3.  The  analysis  presented  in  this  report  predicted  stall- flutter- 
induced  oscillations  of  the  NH-3A  (S-61F)  push-rod  loads  initiating  only 
at  lifts  greater  than  those  at  which  oscillations  were  observed  in  flight. 

4.  Accurate  correlation  of  flatwise  and  edgewise  blade  response  is 
essential,  since  this  response  appears  to  contribute  significantly  to  the 
NH-3A  oscillatory  push-rod  loads,  with  unsteady  aerodynamics  serving  to 
amplify  such  oscillations  in  the  stalled  flow  region.  However,  replacement 
of  steady-state  aerodynamics  with  unsteady  aerodynamics  adversely  affected 
correlation  of  flatwise  and  edgewise  blade  stress  correlation  for  the  NH-3A. 

5.  Analysis  of  the  CH-53A  maneuver  cases  supports  the  NH-3A  con¬ 
clusions  that  both  unsteady  aerodynamics  and  good  flatwise  and  edgewise 
response  predictions  are  needed  to  predict  control  loads  adequately.  In 
addition,  correlation  of  control  loads  during  the  maneuvers  was  more 
encouraging  than  in  steady  flight. 

6.  For  the  flight  conditions  investigated  in  the  Parametric  Studies, 
reducing  the  natural  frequency  of  the  first  torsional  mode  from  8.2P  to  4P 
reduced  the  amplitude  of  the  self-excited  torsional  oscillations  at  a  given 
rotor  lift  because  of  the  slower  response  of  the  4P  blade. 
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7.  Movement  of  the  blade  center -of- gravity  axis  aft  (from  25  percent 
to  jO  percent  chord)  had  i  severe  detrimental  effect  on  blade  torsional 
response  and  associated  oscillatory  aosents.  Uhls  Is  primarily  doe  to  the 
destabilizing  dynamic  coupling  between  the  blade  flapping  and  flatwise 
bending  modes  and  the  torsional  mode.  The  effect  of  aft  center  of  gravity 
on  the  self -excited  moment  components  was  small  but  generally  detrimental. 

8.  Increasing  blade  twist  reduced  the  self -excited  moment  oscilla¬ 
tions  at  a  given  lift  when  rotor  inflow  was  assumed  to  be  constant. 

9-  Except  for  reducing  significantly  the  beneficial  effects  of 
increased  twist  noted  in  (8)  above,  the  inclusion  of  variable  inflow  effects 
had  little  effect  on  the  general  nature  of  the  results  obtained. 

10.  At  conditions  where  retreating  blade  angles  cf  attaO:  exceeded 
steady -state  stall  values,  significantly  higher  rotor  lifts  were  predicted 
using  unsteady  aerodynamic  data  than  were  predicted  using  steady-state 
aerodynamic  data. 
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At  the  present  time,  it  is  difficult  to  establish  the  mason(s )  for 
the  discrepancies  noted  betw  -.a  analytical  and  experimental  results, 
principally  because  of  the  complexity  of  the  experimental  model  used  in  the 
comparison  -  a  full-scale  helicopter.  Achievement  of  good  correlation  in 
such  circumstances  requires  that  (l)  detailed,  accurate  flight  test  data  be 
available,  (2)  all  facets  of  the  aeroelastic  analysis  be  accurate  (inflow, 
applicability  of  sinusoidal  unsteady  data,  use  of  compressibility  scaling,, 
etc.), and  (3)  factors  neglected  in  the  analysis  such  as  radial  flow,  wake 
geometry  distortions,  etc.,  be,  in  fact,  negligible.  Many  approaches  could 
be  followed  in  attempting  to  resolve  the  problem.  One  could,  for 
example,  eliminate  the  need  for  compressibility  scaling  assumptions  by 
acquiring  additional  sinusoidal  unsteady  data  at  higher  Mach  numbers.  It 
is  believed,  however,  that  this  approach  may  be  premature  and  that  the  two 
investigations  recommended  below  represent  the  most  logical  next  steps  in 
the  study  of  stall-flutter-induced  control  loads. 

1.  An  experimental  study  should  be  undertaken,  using  a  two- 
dimensional  airfoil  as  a  model,  to  examine  the  validity  of  an  assumption 
that  forms  the  basis  for  nearly  all  rotor  dynamic  stall  analyses;  namely, 
that  unsteady  aerodynamic  data  obtained  for  an  airfoil  executing  sinusoidal 
motion  can  be  used  to  predict  the  characteristics  of  an  airfoil  executing 
nonsinusoidal  motions  into  stall.  Conclusions  from  such  a  study  would  have 
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far-reaching  effects  on  the  course  of  future  unsteady  aerodynamic  resting 
as  applied  to  helicopter  rotors. 

2.  As  analytical  study  should  be  conducted  to  determine  if  the 
quantitative  accuracy  of  the  method  described  in  this  report  can  be 
Improved  to  permit  its  use  as  an  interim  design  tool  while  more  fundamental, 
longer  range  studies  of  the  type  noted  above  are  pursued.  In  this  study 
consideration  should  be  given  to  possible  seal  empirical  modifications  to 
-ttie  basic  unsteady  data  tabulation  and  associated  compressibility  scaling 
procedures.  In  addition,  the  effect  of  possible  refinements  in  rotor 
trlining  and  inflow  modeling  techniques  on  blade  bending  stress  correlation 
should  be  examined. 
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Figure  1.  Modified  UAC  Rotor  Analysis, 
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NORMAL  FORCE  COEFFICIENT 


Figure  3.  Approximation  of  Typical  Unsteady  Normal  Force  Curve 
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ANGULAR  VELOCITY  PARAMETER,  A 


Figure  k. 


Variation  of  Maximum  Normal  Force  Coefficient,  Lift  Curve  Slope 
and  Angle  for  Zero  Lift  With  Angular  Velocity  and  Angular 
Acceleration  Parameters. 
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Figure  7.  Dynamic  Moment  Characteristic*  -  6  *  +0.01 
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Figure  9*  Comparison  of  Original  and  Reconstructed  Moment  Hysteresis  Loops 
for  Various  Mean  Angles  of  Attack  -  k  =  0.112, 
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Figure  11.  Comparison  of  Original  and  Reconstructed  Moment  Hysteresis 
Loops  for  Various  Reduced  Frequencies  -  aM  =15  deg. 
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Figure  15.  Comparison  of  Original  and  Reconstructed  Normal  Force  K&rs tires  is 
Loops  for  Various  Mean  Angles  of  Attack  -  k  =  0.225. 
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Figure  17.  Steady-State  Normal  Force  Characteristics  for  HACA  0012  Airfoil 
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Figure  18.  Comparison  of  Experimental  Steady- State  Lift  Characteristics 
for  KACA  0012  Airfoil  With  Characteristics  Derived  Using 
Scaling  Procedures. 
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Figure  19*  Comparison  of  Boeing  Experimental  Normal  Force  Hysteresis  Loops 
With  Loops  Derived  Using  Scaling  Procedures  -  M  =  0.4. 
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Figure  20.  Comparison  of  Boeing  Experimental  Normal  Force  Hysteresis  Loops 
With  Loops  Derived  Using  Scaling  Procedures-  M  =  0.6. 
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Figure  21.  Steady-State  Moment  Characteristics  for  NACA  0012  Airfoil 
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Figure  23.  Comparison  of  Boeing  Experimental  Moment  Hysteresis  Loops  With 
Loops  Derived  Using  Scaling  Procedures  -  M  =  0.4. 
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Figure  24.  Comparison  of  Boeing  Experimental  Moment  Hysteresis  Loops  With 
Loops  Derived  Using  Scaling  Procedures  -  M  =0.6. 
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Figure  25 »  Comparison  of  Experimental  Kegativ*  Aerodynamic  Damping  Region 
With  Regions  Derived  Using  Scaling  Procedures. 
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Figure  26.  Comparison  of  Torsional  Mode  Shapes  and  Frequencies  of  Reference 
and  Loir-Frequency  Blades. 
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Figure  27.  Effect  of  Unsteady  Aerodynamics  on  Time  Histories  of  Blade 
Root  Pitching  Moment  -  Reference  Blade. 
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Figure  27.  Concluded. 
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AZIMUTH  ANGLE,  «  ,  DEG 

Figure  28.  Effect  of  Unsteady  Aerodynamics  on  Time  Histories  of  Angle  of 
Attack  at  75  Percent  Blade  Span  -  Reference  Blade. 
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Figure  30.  Effect  of  Variable  Inflow  on  Time  Histories  of  Blade'  Root 
Pitching  Moment  -  Reference  Blade. 
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Figure  30.  Concluded. 
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Figure  U.  Parameters  Defining  Regions  in  Which  Scaled  Unsteady  Aerodynamics  Are  Used 
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Figure  32.  Sensitivity  of  Reference  Blade  Results  to  Mach  Number  Scaling 
Limits,  Variable  Inflow  iterations,  and  Treatment  of  Prandtl- 
Glauert  Factor. 


80 


BLADE  HOOT  PITCHING  MOMENT  COEFFICIENT/SOLIDITY,  (C 


fi=0.4 


- REFERENCE  BLADE  (-V®TWIST) 

- HIGH  TWIST  BLADE  <-*“  TWIST) 


AZIMUTH  ANGLE,  i jj,  DEG 
(a)  $75=  12.2  AND  14.2  DEG 

Figure  33*  Effect  of  Blade  Twist  on  Blade  Root  Pitching  Moment-Constant 
Inflow. 
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Figure  33.  Concluded. 
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Figure  3^ •  Effect  of  Torsional  Frequency  on  Blade  Root  Pitching  Moment  - 
Constant  Inflow. 
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Figure  31*-  Concluded. 
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Figure  35.  Effect  of  Blade  Chordwise  eg  Position  on  Blade  Root  Pitching 
Moment  -  Constant  Inflow. 
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Figure  35*  Concluded. 
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Figure  37.  Effect  of  Blade  Twist  on  Blade  Root  Pitching  Moment  -  Variable 
Inflow. 
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Figure  37*  Concluded 
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Figure  38.  Effect  of  Torsional  Frequency  on  Blade  Root  Pitching  Moment 
Variable  Inflow. 
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Figure  38.  Concluded. 
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Figure  39*  Effect  of  Blade  Chordwise  eg  Position  on  Blade  Root  Pitching 
Moment  -  Variable  Inflow. 
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Figure  4l.  Comparison  of  Blade  Root  Pitching  Moment  and  Flapping  Mode 
Acceleration  -  Constant  Inflow. 
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Figure  4l.  Comparison  of  Blade  Root  Pitching  Moment  and  Flapping  Mode 
Acceleration  -  Constant  Inflow. 
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figure  42.  Comperieoc  of  Experimental  and  Con* tent  Inflow 
Analytical  Puah-Rod  Moment*,  Cut  43. 
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Comparison  of  Experimental,  Variable,  and  Conatant 
Inflow  Analytical  Push-Rod  Moments,  S-61F,  Case  45. 
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Figure  44.  Comparison  of  Experimental  and  Constant  Inflow 
Analytical  Push- Rod  Moments,  S-61F,  Case  42. 
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Figure  45.  Comparison  of  Experimental  and  Constant  Inflow 
Analytical  Push- Rod  Moments,  S-61F,  Case  44. 
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Figure  46.  Comparison  of  Experimental  and  Constant  Inflow 
Analytical  Push-Rod  Moments t  8-61F,  Case  36. 
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Figure  47.  Comparison  of  Experimental  and  Constant  Inflow 
Analytical  Push- Rod  Moments,  S-61F,  Case  35. 
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Figure  48.  Comparison  of  Blade  Flatwise  Air  Loading  Distri¬ 
butions,  S-61F,  Case  43. 
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Figure  49.  Comparison  of  Blade  Flatwise  Air  Loading  Distri¬ 
butions,  S-61F,  Case  45. 
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Figure  50,  Comparison  of  Blade  Flatwise  Air  Loading  Distri¬ 
butions,  S-61F,  Case  42. 
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Figure  51.  Comparison  of  Blade  Flatwise  Air  Loading  Distri 
buttons,  S-61F,  Case  36. 
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Comparison  of  Experimental l  Constant  Inflow,  a 
Modified  Airloads  Analytical  Push- Rod  Moments, 


0009 


Figure  54.  Comparison  of  Experimental ,  Conatant  Inflow,  and 
Modified  Airloads  Analytical  Push-Rod  Momenta, 
3-61F,  Case  42. 


CASE 


Figure  35.  Comparison  of  Experimental ,  Cons  Mot  Inflow,  ami 
Modified  Airload*  Analytical  Fuah-Eod  Momenta, 
S-61F,  Caae  36. 
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Figure  56.  Comparison  of  Experimental  Push- Rod  Moment  With 

Analytical  Push- Rod  Moment  at  Increased  Collective 
Pitch,  S-61F,  Case  43. 
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Figure  57.  Comparison  of  Experimental  Push- Rod  Moment  With 
Product  of  Experimental  Flatwise  and  Edgewise 
Bending  Moments,  S-61F,  Case  43. 
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AZIMUTH  ANGLE, •*»,  DEO 

Figure  58.  Comparison  of  Experimental  Push- Rod  Moment  With 
Analytical  Push-Rod  Moments  for  Steady  and  Un¬ 
steady  Aerodynamic  Data,  S-61F,  Case  43. 
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Figure  62.  Comparison  of  Experimental  and  Analytical  Puah-Rod 
Moments  for  Tenth  Revolution,  CH-53A,  120- Knot, 

60- Degree  Right  Turn. 
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Figure  63.  Time  History  of  Lift,  CH-53A,  150-Knot  Right 
Turn  at  Varying  Angle  of  Bank. 
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Figure  64.  Time  History  of  Push-Sod  Moment  Amplitude, 

CK-53A,  150-Knot  Right  Turn  at  Varying  Angle 
of  Bank. 
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Figure  65.  Comparison  of  Experimental  and  Analytical  Push-Rod 
Moments  for  Twentieth  Revolution,  CH-53A,  150-Knot 
Right  Turn  at  Varying  Angle  of  Bank. 
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Figure  66.  Comparison  of  Experimental  and  Analytical  Push- Rod 

Moments  for  Forty-First  Revolution,  CH-53A,  150-Knot 
Right  Turn  at  Varying  Angle  of  Bank. 
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Figure  69.  Typical  Variation  of  Pitching  Moment  Coefficient  With  Frequency 
for  Mean  Angle  Above  Steady-State  Stall  Angle. 
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Figure  70.  Flow  Chart  for  Data  Conversion  to  a,  A>&nd  B  Form. 
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PITCHING 


ANGULAR  VELOCITY  PARAMETER,  A 

Figure  71 .  Typical  Variation  of  Pitching  Moment  Coefficient  With  Angular 
Velocity  Parameter  and  Comparison  With  Theory. 
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Typical  Variation  of  Normal  Force  Coefficient  With  Angular 
Velocity  Parameter  and  Comparison  With  Theory. 


TABLE  I.  TABULATION  OF  UNSTEADY  NORMAL  FORCE  COEFFICIENT  Cn  FOR  M 
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TABLE  I  -  CONTINUED 


a  =  9  DEGREES 


B 

A=-.04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

•  65 

.65 

■69 

•71 

•  75 

•  85 

.96 

-.004 

.72 

•75 

.79 

.80 

.82 

.89 

•  98 

-.001 

.72 

•75 

.84 

.88 

•90 

•93 

1.03 

0 

.74 

.70 

•87 

•93 

•91 

•95 

1.05 

.001 

•73 

•  69 

.84 

•93 

•91 

•  96 

1.06 

.004 

•71 

•71 

•75 

.74 

•  83 

•96 

1.10 

.01 

.74 

•  70 

.68 

.72 

•78 

•  86 

1.06 

a  a  10  DEGREES 


B 

A=-.04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

•75 

•75 

.80 

.81 

.85 

•  95 

1.05 

-.004 

•  83 

•85 

.88 

.90 

•  93 

1.00 

1.07 

-.001 

.82 

.80 

.86 

■97 

1.01 

1.03 

1.11 

0 

,81 

•76 

.86 

1.02 

1.02 

1.04 

1.14 

.001 

.80 

•71 

•  73 

1.00 

1.00 

1.06 

1.14 

.004 

.74 

•73 

•73 

•  77 

.88 

1.02 

1.18 

.01 

•78 

.68 

•71 

•75 

.81 

.88 

1.11 

a  =  11  DEGREES 


B 

A=-.04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

•85 

.86 

•91 

•93 

•  96 

I.05 

1.14 

-.004 

•93 

•95 

•  98 

1.00 

1.06 

1.10 

1.17 

-.001 

.92 

.86 

•92 

1.05 

1.11 

1.12 

1.20 

0 

•89 

•79 

.34 

1.08 

1.12 

1.12 

1.21 

.001 

.86 

•72 

.80 

1.02 

1.10 

1.14 

1.21 

.004 

.76 

•71 

•  70 

.80 

•93 

1.08 

1.23 

.01 

.80 

•67 

.72 

•79 

.84 

•91 

1.14 
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TABLE  I-  CONTINUED 


a  =  12  DEGREES 


B 

■4- 

o 

i 

ii 

< 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

•95 

•96 

1.01 

1.03 

1.06 

1.14 

1.23 

-.004 

1.03 

1.06 

1.08 

1.10 

1.17 

1.20 

1.27 

-.001 

1.00 

•90 

•  94 

1.12 

1.22 

1.21 

1.27 

0 

•95 

.80 

.80 

1.12 

1.24 

1.21 

i  .29 

.001 

•91 

•72 

•  73 

1.02 

1.19 

1.22 

1.27 

.004 

•77 

.68 

•65 

.81 

•  97 

1.13 

1.26 

.01 

.81 

.67 

•72 

.80 

.86 

•93 

1.17 

a  =  13  DEGREES 


B 

A=- .04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

1.05 

1.06 

1.11 

1.13 

1.16 

1.24 

1-31 

-.004 

1.14 

1-15 

1.18 

1.20 

1.28 

I.30 

1-35 

-.001 

1.08 

•92 

•95 

1.17 

1.32 

I.29 

1.34 

0 

1.00 

.80 

•  76 

1.12 

1-33 

1.28 

1.34 

.001 

•94 

•71 

.66 

1.01 

1.26 

1.28 

1.32 

.004 

•77 

.64 

.61 

.81 

1.01 

1.16 

1.26 

.01 

.82 

.68 

•  73 

.82 

.88 

.96 

1.19 

a  =  14  DEGREES 


B 

A=- . 04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

1.15 

1.17 

1.21 

1.24 

1.26 

1.32 

1.40 

-.004 

1.24 

1.26 

1.28 

1.30 

1.38 

1-39 

1.44 

-.001 

1.14 

•93 

•96 

1.22 

1.41 

1.36 

1.4o 

0 

1.04 

.78 

•74 

1.12 

1.41 

1-35 

1-39 

.001 

.004 

•  96 

.70 

.61 

1.00 

1-32 

1.34 

1-37 

•76 

•59 

•  59 

.81 

1.04 

1.20 

1.26 

.01 

.82 

.68 

•73 

.84 

•  90 

•98 

1.21 

L30 


TABLE  I  -  COffTINUED 


a  =  15  DEGREES 


B 

A— .04 

-.025 

-.01 

0 

.01 

•  .025 

.04 

-.01 

1.26 

1.28 

1-31 

1.34 

1.36 

1.40 

1.48 

-.004 

1-33 

1-35 

1.38 

1.40 

1.48 

1.48 

1-52 

-.001 

1.18 

.92 

•96 

1.25 

1.49 

1.44 

1:46 

0 

1.08 

•76 

•71 

1.11 

1.46 

1.40 

1.43 

.001 

•96 

.68 

.60 

1.00 

1.38 

1-39 

1.40 

.004 

.76 

•57 

•59 

.82 

1.07 

1.22 

1.25 

.01 

.81 

.69 

.74 

.86 

.92 

1.01 

1.23 

a  =  l6  DEGREES 

B 

A=- . 04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

1.36 

1.38 

1.41 

1.44 

1.46 

1.49 

I.56 

-.004 

1.40 

1.43 

1.46 

I.50 

1.58 

1-57 

1.60 

-.001 

1.19 

.91 

•96 

1.28 

1-57 

1.50 

1.51 

0 

1.08 

•75 

.69 

1.10 

1-53 

1.46 

1.46 

.001 

■95 

.65 

.60 

1.00 

1.43 

1.43 

1.43 

.004 

•  75 

.56 

.62 

•  83 

1.09 

1.24 

1.23 

.01 

.80 

.70 

•76 

.88 

.94 

1.04 

1.24 

a  =  17  DEGREES 

B 

> 

11 

1 

b 

-p- 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

1.46 

1.48 

1-51 

1.54 

1.56 

I.58 

1.64 

-.004 

1-39 

1.49 

1-53 

x.58 

1.66 

1.66 

I.67 

-.001 

1-15 

.90 

•96 

1.28 

1.62 

1-57 

1-55 

0 

1.04 

.74 

•70 

1.10 

1.58 

1-53 

I.50 

.001 

•93 

.64 

.62 

1.00 

1.46 

1.47 

1.44 

.004 

.74 

•57 

.65 

86 

1.11 

1.25 

1.20 

.01 

.78 

.71 

•77 

•91 

.96 

1.08 

1-25 

TABLE  I  - 

COHTXNUED 

a  «  18. DEGREES 

B 

A— .04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

l.M> 

1-59 

1.6l 

1.64 

1.66 

1.66 

1-71 

-.004 

1-35 

1.47 

1.58 

1.64 

1.74 

1.74 

1.74 

-.001 

1.08 

.89 

•95 

1-27 

1.66 

1.64 

1.60 

0  * 

1.00 

•75 

•70 

1.10 

1.60 

1.60 

1.54 

.001 

•90 

.64 

.65 

1.00 

1.47 

1.54 

1.47 

.004 

•73 

.60 

.68 

.88 

1.13 

I.27 

1.22 

.01 

76 

.72 

•79 

■93 

•99 

1.12 

1.26 

a  -  20  DEGREES 

B 

A— .04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

1.77 

1-79 

1.80 

1.84 

1.86 

1.84 

1.87 

-.004 

1.17 

1.25 

1.56 

1-71 

1.87 

1-93 

1.91 

-.001 

•90 

.84 

■93 

1.21 

1.70 

1.82 

1*75 

0 

.84 

.73 

•75 

1.08 

1.61 

1-73 

1.67 

.001 

.83 

.66 

•  71 

1.00 

1.47 

1.64 

1-59 

.004 

•72 

.68 

•  76 

•92 

1.16 

1.33 

1.41 

.01 

•72 

.74 

• 

.84 

•99 

1.06 

1.21 

1-31 

1  *  12  DEGREES 

-.04 

-.025 

..c> 

0 

.01 

.025 

.04 

L.91 

1-95 

1-99 

2.02 

2.06 

2.05 

2.07 

L.07 

1.08 

1.47 

1.70 

1.95 

2.13 

2.12 

.82 

.78 

•92 

1.16 

1.71 

I.98 

I.98 

•78 

.72 

•79 

1.07 

1.53 

1.85 

1.91 

•76 

.70 

•79 

1.02 

1.42 

1.72 

1.84 

.72 

•75 

.85 

•99 

1.18 

1.38 

1.61 

.72 

.81 

•  92 

1.04 

1.14 

1-31 

1.4i 

table  I- concluded 


a  =  24  DEGREES 


B 

A=-.04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

1.85 

1.99 

2.16 

2.20 

2.23 

2.28 

2.25 

-.004 

1.05 

i.05' 

'  1.35 

1.60 

1.96 

2.27 

2.27 

-.001 

.86 

.80 

•91 

1.15 

1.58 

1-95 

2.10 

0 

■83 

.78 

.84 

1.08 

1.44 

1.84 

2.04 

.001 

.78 

.76 

•85 

1.05 

1-35 

1.69 

1.95 

.004 

.74 

.81 

.94 

1.04 

1.20 

1.40 

1.68 

.01 

•76 

•93 

1.05 

1.09 

1.25 

1-37 

1.50 

a  -  26  DEGREES 

B 

A—. 04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

1.67 

1.00 

2.02 

2.14 

2.21 

2.34 

2.35 

-.004 

1.06 

1.07 

1.22 

1.42 

1.83 

2.17 

2.32 

-.001 

•  90 

•87 

.94 

1.19 

1.48 

1-79 

2.01 

0 

.86 

•85 

•  90 

1.14 

1.40 

1.64 

1.84 

.001 

.82 

.82 

■91 

1.11 

1-34 

1.54 

1.68 

.004 

•  76 

•  85 

•99 

1.06 

1.20 

1-35 

1-52 

.01 

.84 

•96 

1.02 

1.07 

1.22 

1-33 

1-56 

a  =  27  DEGREES 

B 

A— .04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

1-55 

1.58 

1.69 

1.88 

2.07 

2.30 

2.38 

-.004 

1.06 

1.09 

1.18 

1.36 

1.70 

2.03 

2.31 

-.001 

.92 

•91 

•98 

1.22 

1.47 

1.72 

1.91 

0 

.86 

.88 

.94 

1.18 

1.42 

1-57 

1.75 

.001 

.85 

•85 

.94 

1.14 

1.36 

1.49 

1.63 

.004 

•79 

•85 

•96 

1.03 

1.19 

I.29 

1.46 

.01 

•  89 

•  90 

•95 

1.04 

1.14 

I.27 

1.59 
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TABLE  II.  TABULATION  OF  UNSTEADY  PITCHING  MCMENT 
COEFFICIENT  Cm(;/4  FOR  M  =  0 


a  =  0  THROUGH  6  DBGREES 

B 

A=-.04 

-.025 

-.01 

0 

J'l 

.025 

.04 

-.01 

.069 

.045 

.022 

.006 

-.010 

-.033 

-.057 

-.007 

.067 

.043 

.020 

.004 

-.012 

-.035 

-.059 

-.004 

.065 

.042 

.018 

.002 

-,013 

-.037 

-.061 

-.001 

.063 

.040 

.016 

.001 

-.015 

-•039 

-.062 

0 

.063 

•039 

.016 

0.0 

-.016 

-.039 

-.063 

.001 

.062 

•039 

.015 

-.001 

-.016 

-.040 

-.063 

.00 4 

.061 

.037 

.013 

-.002 

-.018 

-.042 

-.065 

.007 

.059 

•035 

.012 

-.004 

-.020 

-.043 

-.067 

.01 

.057 

•033 

.010 

-.006 

-.022 

-.045 

-.069 

a  = 

7  DEGREES 

B 

A=-.04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

.066 

.043 

,020 

.003 

-.013 

-.035 

-.058 

-.007 

.064 

.040 

.016 

.001 

-.013 

-.033 

-.062 

-.004 

.057 

•035 

.015 

0.0 

-.015 

-.039 

*.065 

-.001 

.056 

■033 

.010 

-.002 

-.018 

-.043 

-.065 

0 

.050 

.031 

.011 

0.0 

-.015 

-.041 

-.067 

.001 

.050 

.030 

.010 

-.003 

-.018 

-.046 

-.066 

.004 

-035 

.020 

.010 

-.004 

-.020 

-.044 

-.067 

.007 

.016 

-013 

.007 

-.008 

-.023 

-.051 

-.074 

.01 

-.010 

-.002 

0. 

-.003 

-.026 

-.050 

-.076 

a  = 

8  DEGREES 

B 

A=-.04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

•059 

.040 

• 

0 

►— 1 

CT\ 

0.0 

-.015 

-.040 

-.061 

-.007 

.056 

.034 

.010 

-.002 

-.016 

-.043 

-.065 

-.004 

.045 

.025 

.010 

-.005 

-.017 

-.043 

-.067 

-.001 

.045 

.024 

0.0 

-.005 

-.021 

-.045 

-.068 

3 

.031 

.020 

0.0 

0.0 

-.016 

-.045 

-.071 

.001 

.028 

.010 

0.0 

-.007 

-.024 

-.050 

-.070 

.004 

-.010 

-.016 

.007 

-.006 

-.020 

-.056 

-.070 

.007 

-.015 

-.007 

0.0 

-.010 

-.023 

-.054 

-.074 

.01 

-.030 

-.025 

-.010 

-.012 

-.028 

-.052 

-.080 

TABLE  II  -  CONTINUED 

a  = 

9  DEGREES 

B 

A=-.04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

.043 

•035 

.010 

-.003 

CO 

r— 1 

0 

• 

1 

-.045 

-.067 

-.007 

.044 

.026 

.006 

-.007 

-.020 

-.047 

-.071 

-.004 

•  038 

.020 

.005 

-.011 

-.021 

-.048 

-.073 

-.001 

.035 

.015 

-.005 

-.012 

-.024 

-.050 

-.074 

0 

.020 

.007 

-.005 

0.0 

-.016 

-.049 

-.076 

.001 

.013 

-.025 

-.005 

-.012 

-.034 

-.055. 

-.074 

.004 

-.028 

-.027 

0.0 

-.009 

-.020 

-.049 

-.074 

.007 

-.034 

-.030 

-.007 

-.011 

-d 

CM 

O 

• 

1 

-.054 

-.074 

.01 

-.0>4 

-.044 

-.021 

-.018 

-.029 

-.053 

-.080 

a  - 

10  DEGREES 

B 

A--.04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

•  039 

.029 

.004 

-.007 

-.024 

-.049 

-.073 

-.007 

.033 

-023 

.002 

-.011 

-.024 

-.050 

-.075 

-.004 

-035 

.016 

0.0 

-.016 

-.027 

-.053 

-.077 

-.001 

.025 

.010 

-.009 

-.017 

-.027 

-.054 

-.080 

0 

.005 

-.004 

-.015 

0.0 

-.017 

-.049 

-.081 

.001 

.016 

-.040 

0.0 

-.020 

-.04l 

-.057 

-.081 

.004 

-.047 

-.035 

-.002 

-.010 

-.020 

-.050 

-.076 

.007 

-.054 

-.045 

-.015 

-.010 

-.025 

-.053 

-.075 

.01 

-.075 

-.060 

-•033 

-.025 

-.028 

-.052 

-.078 

a  = 

11  DEGREES 

P 

A=-.04 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

•033 

.023 

-.001 

-.012 

-.027 

-.053 

-.077 

-.007 

.036 

.021 

.001 

-.014 

-.027 

-.053 

-.077 

-.004 

.032 

.015 

-.003 

-.018 

-.032 

-.055 

-.080 

-.001 

.010 

.002 

-.012 

-.020 

-.030 

-.056 

-.084 

0 

-.020 

-.014 

-.025 

0.0 

-.020 

-.054 

-.087 

.001 

-.005 

-.035 

0.0 

-.023 

-.044 

-.071 

-  .086 

.004 

-.070 

-.040 

-.005 

-.010 

-.020 

-.050 

-.078 

.007 

-.078 

-.053 

-.025 

-.022 

0.025 

-.052 

-.075 

.01 

-.082 

-.072 

-.050 

-.035 

-.027 

-.048 

-.078 
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a  = 

12  DEGREES 

B 

A=-.04 

IT\ 

CVI 

O 

• 

1 

-.01 

0 

.01 

.025 

.04 

-.01 

.031 

.016 

-.003 

-.015 

-.030 

-.054 

-.080 

-.007 

.035 

.020 

0.0 

-.015 

-.030 

-.055 

-.080 

-.00!* 

.030 

.01? 

-.005 

-.020 

-.035 

-.057 

-.083 

-.001 

-.008 

-.010 

-.013 

-.022 

-.031 

-.058 

-.085 

0 

-.01*0 

-.028 

-.033 

0.0 

-.023 

-.057 

-.089 

.001 

-.011 

-.014 

.007 

-.019 

-.043 

-.087 

-.088 

.001* 

-.090 

-.040 

-.005 

-.007 

-.020 

-.048 

-.079 

.00? 

-.096 

-.062 

-.028 

-.025 

-.023 

-.038 

-.072 

.01 

-.100 

-.0S0 

-.050 

-.04c 

-.025 

-.040 

-.065 

a  = 

13  DEGREES 

B 

A=-.04 

-.025 

-.01 

0 

.01 

.025 

.04 

a 

• 

o 

.031* 

.019 

0.0 

-.015 

-.030 

-.055 

-.080 

-.007 

.031* 

.020 

0.0 

-.016 

-.030 

-.055 

-.082 

-.001* 

.011 

.014 

-.007 

-.020 

-.038 

-.060 

-.085 

-.001 

-.036 

-.046 

-.018 

-.023 

-.030 

-.060 

-.085 

0 

-.050 

-.050 

-.041 

0.0 

-.028 

-.OoO 

-.067 

.001 

-.050 

-.031 

.010 

0.0 

-.035 

-.084 

-.084 

.001* 

-.100 

-.042 

-.005 

0.0 

-.015 

-.045 

-.079 

.007 

-.108 

-.070 

-.027 

-.010 

1 

• 

O 

vn 

-.033 

-.060 

.01 

-.108 

-.072 

-.o4o 

-.021 

-.020 

-.037 

-.057 

a  = 

14  DEGREES 

B 

A=-.0l* 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

.032 

.022 

.003 

-.016 

-.031 

-.055 

-.083 

-.007 

.025 

.020 

0.0 

-.016 

-.032 

-.057 

-.083 

-.004 

.030 

0.0 

-.010 

-.021 

-.038 

-.060 

-.087 

-.001 

-.080 

-•095 

-.032 

-.026 

-.032 

-.056 

-.085 

0 

-.090 

-.072 

-.047 

-.005 

-.036 

-.060 

-.083 

.001 

-.104 

-.046 

.010 

.006 

-.028 

-.082 

-.080 

.001* 

-.112 

-.046 

-.007 

0.0 

-.012 

-.042 

-.076 

.007 

-.119 

-.070 

-.026 

-.014 

-.017 

-.032 

-.057 

.01 

-.106 

-.057 

-.036 

-.020 

-.020 

-.023 

-.050 

TABLE  n  -  CCROTJaD 


a  = 

17  DEGREES 

B 

3 

• 

1 

ii 

-.025 

-.01 

0 

.01 

.025 

.04 

-.01 

-.030 

-.011 

-.010 

-.020 

-.033 

-.060 

-.085 

-.007 

-.105 

-.065 

-.034 

-.020 

-.036 

-.069 

-.093 

-.004 

-.145 

-.120 

-.055 

-.020 

-.034 

-.061 

-.100 

-.001 

-.164 

-.135 

-.068 

-.041 

-.034 

-.074 

-.106 

0 

-.157 

-.09  6 

-.054 

-.080 

-.057 

-.051 

-.063 

.001 

-.145 

-.056 

-.007 

-.025 

-.016 

-.043 

-.054 

.004 

-.126 

-.076 

-.055 

-.055 

-.060 

-.048 

-.039 

.00? 

-.090 

-.091 

-.097 

-.106 

-.100 

-.074 

-.031 

.01 

-.079 

-  .068 

-.070 

-.061 

-.075 

-.065 

-.030 
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TABLE  III. 

SCALING  PARAMETER  VALUES  FOR  BOEING  AIRFOIL 

Mach  Number 

Normal  Force 
Stall  Angle 

Normal  Force 

Coef.  at  Stall 

Pitching  Moment 
Stall  Angle 

M 

r 

n  STALL ,ss 

“sm>  de6 

0.4 

8.4 

0.85 

10.0 

0.6 


4.8 


0.53 


7.0 


.ABLE  IV.  SCALING  PARAMETER  VALUES  FOR  UAC  AIRFOIL 


tech  N  mber 

M 

Normal  Farce 
Stall  Angle 
a$n»  dee 

Normal  Force 

Coef.  at  Stall 

CnSTALL,SS 

Original 

Pitching 

Moment 

Stall  Angle 
a*m>  deg 

Revised 

Pitching 

Moment 

Stall  Angle 
a*m  ,  deg 

0 

10.0 

1.05 

13.5 

13-5 

0.2 

10.0 

1.05 

13.5 

13.5 

0.3 

10.4 

1.09 

13.0 

12.5 

0.4 

8.4 

0.85 

12.3 

11.0 

0.5 

6.S 

0.69 

11.5 

9.0 

0.6 

4.8 

0.53 

10.6 

7.0 

0.7 

2.9 

0.37 

10.6 

7.0 

0.8 

2.0 

0.23 

10.6 

7.0 

0.85 

2.0 

0.23 

10.6 

7.0 
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TABLE  V. 

BLADE  DESIGNS 

i _  Z2 

Design 

Twist  Rate, 
B  |  ,  deg 

Torsional 
Frequency „ 

\ 

Chordwise  eg 
Location, 

y'Oc«/C 

1. 

Reference  Blade 

-4 

8.2P 

0.25 

2. 

High-Twist  Blade 

-8 

8.2P 

0.25 

3- 

Low-Frequency  Blade 

-4 

4p 

0.25 

4. 

Afty-cg  Blade 

-4 

8.2P 

0.30 
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TABLE  VI.  REFERENCE  31/DE  CHARACTERISTICS 


I. 

Chord  (ft^ 

1.52 

2. 

Twist  rate  (deg) 

-4 

3- 

Radius  (ft) 

31 

4. 

Flap  -  lag  hinge  offset 

0.054 

5. 

Tip  speed  { ft/sec) 

662 

6. 

Pitch-flap  coupling  ratio 

0.17 

T  • 

Flatwise  bending  frequency  ratios 

2.66P,  4.87P,  7.68P 

8. 

Edgewise  bending  frequency  ratios 

3.44p,  9-OP 

9. 

Torsional  frequency  ratios 

S.22P,  21.96P 

10. 

Control  system  stiffness 
(ft-lb/rad) 

1.2  x  106 

11. 

Nominal  section  weight  (lb/in.) 

0.5 

12. 

Nominal  section  flatwise  stiffness 
(lb/ln.2) 

2.2  x  107 

13. 

Nominal  section  edgewise  stiffness 
(lb/in.2) 

24  x  107 

14. 

Nominal  section  torsional  stiffness 
(lb- in.2) 

2.9  x  107 

15. 

Nominal  section  torsional  Inertia 
(lb-se^2) 

0.024 

16. 

Section  eg  (percent  chord) 

25 

TABLE  VII 

.  SUMMARY  OF 

S-61F  FLIGHT  TEST  CASES  USED  FOR  OOHRELATIC* 

1  ■ 

Case 

Velocity, 

Rotor 

Boaioal 

logical  Throat 

Staber 

Knots 

Axial  Force, 

Advance 

Coefficient- 

Pounds 

Ratio 

Solidity  Ratio 

35 

160.80 

5060 

.*10 

.020c 

36 

166.50 

13300 

•  *25 

.05*3 

*2 

l86.*5 

8630 

.*?6 

.0351 

*3 

185.97 

1*200 

.0579 

kli 

184.07 

*290 

.*69 

.0175 

*5 

186.92 

9870 

.*77 

.0*02 
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APPENDIX  I 

MCDIFICA  TICKS  TO  EglATIOKS  CF  BHTSHKE  13 


This  appendix  is  intended  to  be  used  in  conjunction  with  References 
13  and  lit  and,  when  so  used,  to  provide  documentation  of  the  specific 
modifications  made  under  this  contract  to  the  blfde  equations  of  notions 
given  in  those  references .  Inasmuch  as  Reference  13  and  14  are  generally 
available  and  would  be  required  for  a  detailed  understanding  of  the  blade 
equations  used  in  this  study,  no  attempt  has  been  Bade  to  reproduce  herein 
any  parts  of  References  13  and  14.  As  a  result,  cany  of  the  symbols  In 
the  following  equations  are  net  Included  In  the  List  of  Symbols  of  this 
report  but,  rather,  appear  in  the  List  of  Symbols  of  References  13  and  14. 

Tile  basic  equations  governing  the  motion  of  rotor  blades  having  chord- 
vise  mass  unbalance  and  operating  in  steady,  level  flight  are  given  in 
Reference  13-  To  extend  these  to  the  case  of  maneuvering  flight,  the 
equations  were  modified  by  adding  terms  representing  the  effect  of  fuselage 
motion  on  blade  dynamic  and  aerodynamic  forces .  These  terms  were  obtained 
directly  from  Reference  14  and,  for  convenience,  are  given  below  for  each 
of  the  equations  of  Reference  13  to  which  they  apply. 

Flatwise  Equation  of  Motion 

Ob  the  right-hand  side  of  Eq.  (80)  of  Reference  13  is  added 
C,3j(Asine0+ «Z(COS0O)  +  (c,2j-0.75C,0.)  0,  (csin0o  +  Bcos  0O) 

(21) 

+  C,2j  (-  G  cos 60  +  F"sin  9C)  +  C|0.  Ecos0o  +  D  sin£?0) 


Edgewise  Equation  of  Motion 

To  the  right-hand  side  of  Eq.  (86)  of  Reference  13  is  added 

l  jt_ 

C33p(-Acos0o  +  sin  90j  +  (c47p-O.75C,5p)0,  (gsin 0O  -  C cos0o 
“  C|5p  (dcos 9q  +  Esine0)  -  c47p  (f cos 60  +Gsin0o) 
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(e  + r")(wy  cos  f-  wxsinf|  +  £v0y  (sin*  ♦  Boost) 

-78  (w^cosf  +  5y(sin^f|  +  v0^  |sint  +  Scost)  jw#'cos 8  +  ve'sin®|(27) 
+  0,  v0  sin  f  (ve  cos  6  -  «g  sinfi) 


In  Eqs.  (21)  throjgh  (25),  the  following  definitions  apply: 


a  =  wS|sint  -  oncost  +  2(s,,(C0St  +  C,(sint) 


(28) 


B  =  ^  cost  -  ^sint  -  cost  +  «**,,sin  t) 
+  wI(  <7^  sint  +  vQx  cost) 


(29) 


C  s  %,t  '  \wh  +VOyiW», 


(30) 


D  =  -r 


(3D 


r  cost  -  -r  sin  t  -  v  (w  cost  +  wy.sint) 

*  (70y  sin  t  +  v0i  cos  t )  +  e  ^  -  8  ( ^  cos  t  +  -fc  sin  t ) 

1  1  *i  °yi 

-  "vo,  8  (w  cost  -  <2,  sint)  +  o.  8(-v0  sint+70  cost) 

I  I  I|  Jfi 
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E  *  \  v0  ♦  Bk,  v0  -  V  0  cos  *  *  sin*> 

*i  •  >i  *i 

+  21(3*  sin  *  +  5^  cos  *  )  -  e  (*3^  cos  *  - -^  sin*) 
-/3(-r  cos*+-$r  sin*)  +  £5.  cos*  -  70  sin*) 


(32) 


F  =  *5*ti  -  2/3  («,,  cos  *  +  s|(sin*)  -  0{  cos  *  +  sin  * )  (33) 


Z  *  sin*  -  "ST  cos  *  +  2(3.  cos  *  +  5,  sin  *)  +  205. 

•  #1  I  1  ' 

+  28  (ot|co*  *  ♦  5S,(  sin  *)  +  8  (  3^  cos  *  +  sin  *)  (3*0 

+  28  (o,(  cos  *  -  5k, sin  *  ) 


In  addition  to  these  modifications  to  the  equations  of  Reference  13, 
other  changes  to  the  basic  aerodynamic  equations  were  necessary  when 
unsteady  data  were  used.  Specifically,  when  a  blade  section  was  operating 
at  a  combination  of  a  Mach  number  and  an  angle  of  attack  where  the  unsteady 
data  were  to  be  used,  the  section  angle  of  attack  was  defined  by  the  veloci¬ 
ty  components  at  the  quarter- chord  rather  than  those  at  the  three-quarter 
chord  point,  and  the  damping  moment  In  pitch  predicted  by  quasi-steady 
theory  was  eliminated.  Formally,  these  changes  were  accomplished  by 
setting  Yioj/kc  e1ual  zero  ln  Eqs .  (139)  and  (lto)  of  Reference  13  and 
Oq  equal  to  0.5  in  Eq.  (1^6). 
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APPHTDU  II 

I>EEAII£  OF  THE  COWVERSIOB  CF  SINUSOIDAL  DATA  TO  a,  A  ,  AND  B  FOBI 


In  any  tabulation  of  data  for  subsequent  use  in  a  digital  computer,  It 
is  imperative  that  the  dependent  quantity  be  associated  vith  convenient 
tabular  values  of  the  Independent  parameters.  In  the  present  case,  the 
dependent  quantity  was  either  cn  or  cme/t  and  the  independent  parameters 
were  a  ,  a  ,  and  B .  In  this  appendix,  the  method  used  for  establishing  this 
dependency  of  Cm^,  on  a,  A,  and  B  will  be  described  in  detail. 
(Substantially  the  same  method  was  used  in  Reference  10  to  relate  Cn  too, 
A,  and  B,  but  this  earlier  report  lacks  the  detail  of  the  present  work.) 

First,  it  is  convenient  to  rewrite  the  reduced  frequency  of  Eq.  (14) 
as 


k  •- 


C  M  .  wcf 

2U  "  U 


(35) 


and  to  inseit  this  form  into  Eqs.  (12)  and  (13)  for  A  and  E  .  Die  result  Is 


A  *  ± 


(36) 


B=  -(^)%-««») 


(37) 


Now  consider  a  typical  hysteresis  loop  of  Cmc/4  versus  a,  as  shown  in 
the  sketch  in  Figure  67.  Die  angular  velocity  parameter  A  is  positive  over 
the  lower  half  of  the  loop,  which  is  the  region  in  which  o  is  increasing 
with  time.  A  is  negative  over  the  upper  half  of  the  loop.  A  reaches  its 
positive  and  negative  maximums  at  a  =  aM  and  vanishes  at  both  angular 
extremes  of  the  loop.  Die  angular  acceleration  parameter  B  is  in  exact 
antiphase  with  the  angular  displacement  from  the  mean  (  a  -  aM),  reaching 
its  negative  maximum  at  maximum  a  and  reaching  its  positive  maximum  at 
minimum  a.  B  vanishes  whenever  a  =  aM.  It  can  be  shown,  both  from  Eqs. 
(36)  and  (37)  and  from  a  consideration  of  the  typical  loop  in  Figure  67, 


* 
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that  convenient  tabular  values  of  a  ,  A ,  and  B  cannot  be  obtained  directly 
from  the  loops.  Hence,  it  is  necessary  to  manipulate  both  the  equations  and 
the  curves  to  produce  cross  plots  suitable  for  interpolation.  In 
subsequent  paragraphs,  each  parameter  in  Eqs.  (36)  and  (37)  will  be  examined 
to  determine  its  suitability  as  an  interpolation  variable. 

The  free  parameters  in  A  and  B  (Eqs.  (36)  and  (37))  are  the  chord  C  , 
the  amplitude  of  motion  a  ,  the  velocity  U  ,  the  angle  of  attack  a  ,  the 
mean  angle  of  attack  aM,  and  the  frequency  f  .  (The  local  angle  of  attack 
a  is  not  a  free  parameter  since  it  has  been  selected  as  one  of  the  three 
tabular  quantities.)  The  chord  C  may  be  eliminated  immediately  because 
only  a  single  value  was  available  from  the  experiment,  C  =  2  ft.  There 
were  only  three  values  of  the  amplitude  o’  and,  as  discussed  in  the  main  body 
of  this  report,  most  of  the  data  were  taken  at  a  =  ±6  deg  and  *8  deg. 
Similarly,  there  were  only  three  values  of  the  velocity  U  (or  Mach  number), 
and  most  of  the  data  were  taken  at  the  lowest  two  values.  The  mean  angle 
of  attack  aM  was  taken  every  three  degrees  over  a  wide  range  of  values . 
However,  two  factors  must  be  considered  here:  (l)  only  integral  values  of 
a  may  be  used  as  mesh  points,  and  (2)  the  displacement  relative  to  the 
mean  position  (a  -aM)  may  not  exceed  the  motion  amplitude  a  .  Hence, 
this  limits  the  number  of  values  of  au  that  could  be  used. 

Ihe  only  remaining  parameter  in  Eqs.  (36)  and  (37)  is  the  frequency  f . 
In  most  cases,  the  frequency  values  at  which  data  were  taken  were  f  =  0 
(steady  state),  4,  8,  12,  and  16  cps  .  This  is  a  set  of  five  equally  spaced 
points  that  extend  over  a  physically  useful  and  realistic  range.  In  many 
cases  in  the  vicinity  of  the  steady-state  stall  angle,  the  three  additional 
values  of  t  =  0.5,  1,  and  2  cps  were  taken.  Hence,  it  is  logical  to  use  f 
as  an  interpolation  variable.  This  prospect  is  made  even  more  attractive 
by  the  manner  in  which  f  appears  explicitly  in  Eqs.  (36)  and  (37)* 

To  demonstrate  the  method  used  herein,  first  solve  Eq.  (37)  for  f  to 
obtain 


f 


(38) 


Wnen  Eq.  (38)  is  substituted  into  Eq.  (38),  the  equation  for  A  becomes 


(39) 


Hie  first  step  in  reducing  the  data  was  to  cross-plot  Cmc>4  versus  f  . 
The  procedure  followed  here  was  to  pick  a  value  of  each  of  the  parameters 
U  (or  M),  a  ,  and  aM  ,  and  then  to  work  with  the  set  of  loops  associated 
with  the  values  of  f  within  the  set.  Two  cross  plots  were  made  through  the 
set  for  each  integral  value  of  a  :  one  cross  plot  was  made  for  the  increasing 
a  branch  of  the  loop,  and  the  other  was  made  for  the  decreasing  a  branch 
of  the  loop.  A  few  selected  cross  plots  are  shown  in  Figure  68  for  a*  « 

0  deg  (representing  the  potential  flow  regime)  and  in  Figure  69  for  = 

15  deg  (representing  the  stalled  flow  regime).  In  both  figures,  M  *  0.2 
and  5  =  6  deg.  It  is  intei-esting  to  note  that  in  potential  flew  (Figure 
68)  the  Cmc/4  variation  with  f  is  well  behaved  over  the  entire  range,  which 
justifies  a  minimal  number  of  frequency  points.  However,  in  stalled  flow 
(Figure  69)  the  Cmc/4  variation  with  f  is  very  erratic,  particularly  in 
the  range  0  <  f  £  b  cps .  These  facts  should  be  noted  in  the  design  of  any 
future  experiments  on  the  unsteady  response  of  oscillating  airfoils . 

After  the  cross  plots  were  completed  for  all  of  the  loops,  the  result 
was  a  multitude  of  curves  of  Cmc/4  versus  f  .  In  this  set  o*-  curves,  all 
values  of  the  parameters  a  ,  aM,  a,  and  U  were  represented  within  the 
limitations  of  the  experiment.  Each  individual  curve  was  associated  with 
a  particular  set  of  values  of  a  ,  aM ,  a  ,  and  U,  and  with  etcher  d>  0  or 
a  <  0. 


In  the  next  part  of  the  data  reduction  procedure,  a  number  of  inter¬ 
related  steps  vere  taker  (see  the  flow  chart  in  Figure  TO) .  Initially,  a 
value  of  a  and  a  value  of  U  were  chosen  (e.g.,  c  =  6  deg,  U  =  230  ft/sec). 
Then  an  integral  mesh  point  value  of  a  was  chosen.  Next,  a  convenient 
mesh  point  value  of  B  was  chosen.  (The  set  of  B  values  used  in  the  present 
study  was  B  =  0,  -0.001,  ±0.004,  ^0 .007,  -0.01.)  At  this  point,  the  only 
undetermined  free  parameter  in  Eqs .  (38)  and  (39),  or  in  the  family  of 
curves  of  Cmc/4  versus  f  ,  is  the  mean  angle  of  attack  aM.  A  series  of 
values  of  aM  was  chosen  next,  subject  to  the  restriction  that  |a  —  aM|< "5  . 
For  each  value  of  aM,  a  value  of  f  was  calculated  from  Eq.  (38)  and  a  value 
of  A  was  calculated  from  Eq.  (39)-  Next,  the  two  curves  of  Crr*,*  versus  f 
for  the  specific  values  of  "a  ,  U,  a,  and  aM  were  selected  from  the  set,  one 
for  increasing  a  and  one  for  decreasing  a .  These  curves  were  entered  at 
the  calculated  value  of  f  ,  and  the  interpolated  values  of  Cm  veT*e 
plotted  versus  the  calculated  value  of  A  .  The  value  of  Cmc/*  associated 


» 
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with  increasing  fl  was  plotted  et  the  positive  value  of  A,  and  the  value  of 
cmc/4  associated  with  decreasing  a  was  platted  at  the  negative  value  of  A . 
For  the  special  case  of  a  stationary  (i.e.,  at  the  end  points  of  a  loop), 
the  interpolated  Cmw  was  plotted  at  A  =  o.  After  this  had  been  done  for 
all  possible  values  of  aM,  a  faired  curve  was  passed  through  the  points, 
representing  the  variation  of  Cme/4  with  A  for  a  single  value  of  a  and  B  • 
A  new  value  of  B  was  then  souglrt,  and  the  entire  process  was  repeated. 


An  example  of  these  curves  is  shown  in  Figure  71  for  M  =  0.2  an--!  fcr 
a  mesh  point  value  of  a  =  IS  deg.  In  this  figure  a  set  of  values  of  B 
ranging  from  -0.007  to  -+0.007  was  selected.  As  discussed  in  the  text,  mcot 
of  the  data  were  originally  obtained  at  5  =  *8  deg;  hence,  the  faired 
curves  are  biased  to  conform  to  these  points .  The  theoretical  variation  of 
Cmc/4  versus  A ,  which  is  derived  in  Appendix  m,  is  also  plotted  in  this 
figure.  It  is  interesting  to  note  that  even  at  this  high  incidence  angle, 
sane  measure  of  agreement  still  exists  between  theory  and  experiment  for 
large  positive  values  of  A  (i.e.,  for  a>  0).  A  study  of  any  of  the  actual 
loops  described  elsewhere  In  this  report  will  show  that  almost  the  entire 
branch  of  the  loop  far  a  >  0  tends  to  behave  like  a  theoretical  loop,  even 
when  the  local  angle  exceeds  the  stalling  angle,  provided  that  the 
frequency  of  motion  Is  high  enough. 

The  faired  curves  of  Cm£i#  versus  A  were  then  cross-plotted  with 
•  respect  to  B  at  constant  a  for  a  selected  set  of  values  of  A ,  namely,  A  = 

0,  id .01,  io.025,  to .04.  Once  again  curves  were  faired  through  the  points, 
and  adjustments  were  made  in  the  faired  curves  to  provide  as  smooth  a  set  as 
possible.  A  third  and  final  cross  plot  was  then  made  of  Cmc/«  versus  a 
for  constant  values  of  A  and  B  (see  Figures  5  through  7)  •  After  the  curves 
were  smoothed  through  these  points,  a  number  of  correlations  were  performed 
between  the  original  loops  and  the  faired  curves.  In  many  cases,  a  normal 
scatter  was  observed  in  correlating  the  smoothed  data  with  the  original 
loops,  but  in  a  few  cases  a  consistent  discrepancy  was  observed.  In  these 
cases,  the  faired  data  points  and  curves  were  adjusted  to  produce  better 
agreement  with  the  original  loops .  This  level  of  agreement  is  discussed  in 
the  main  body  of  this  report  (Figures  9  through  14).  Finally,  after  a 
satisfactory  agreement  had  been  attained,  the  data  points  were  tabulated 
and  are  presented  in  7&ble  II. 
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APP3TDIX  m 

THEORETICAL  K0R4AL  FORCE  AKD  MCXHJT 
COEFFICianS  FOR  SPiUSOIDAI  MOTIOW  IS  TERC  CF  a,  A ,  ASP  B 


The  expressions  for  the  complex,  unsteady  normal  force  and  moment 
coefficients  for  sinusoidal  pitching  notion  will  be  taken  from  Eqs.  (fc) 
and  (5).  The  Theodors en  function,  c(k),  will  be  substituted  from  Eq.  (6), 
and  the  complex  angular  displacement  will  be  taken  from  2q.  (7).  In 
accordance  with  the  statements  in  the  text  following  £q.  (j),  the  first 
terms  of  Eqs .  (k)  and  (5)  will  be  rewritten  as 


2  t  C  (hie  ■  -—2  r  'ly  +  2  ir  C  (k)a  «'"* 


(*» 


*  (°0+i) Cikia* — — » («0+i)a»* ♦  w(°0+?)  (m) 


Hence,  Eqs.  (U)  and  (5)  become 


c.*  =  2  +  {  2(F+,Gi  +  -gfr-  [i+2{F«HG)(^ -o0)]+«fe(|g)2}  (k2) 
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In  the  present  study  the  pivot  axis  was  located  at  the  quarter-chord,  or 
Q0=  -  l/2.  When  this  value  is  substituted  into  these  equations,  the 
complex  exponential  is  expanded,  and  the  real  part  of  each  equation  is 
taken;  the  result  is 
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where 
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is  the  reduced  frequency  for  sinusoidal  motion,  From  Eqs.  (8)  through  (13) 
it  is  easily  shown  that,  for  sinusoidal  motion. 
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Hence,  Bqs .  (4k)  and  (V>)  tay  be  reduced  to  the  fora 
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Although  the  torrional  stall  flutter  phenomenon  is  primarily  associ¬ 
ated  with  the  behavior  of  the  unsteady  pitching  moment,  it  is  instructive 
to  look  at  the  variations  in  normal  force  coefficient  with  changes  in  the 
relevant  parameters.  One  such  variation  is  shown  in  Figure  72,  in  which Cn 
is  plotted  versus  A  for  the  special  case  of  B  =  0.  In  calculating  this 
theoretical  curve,  the  angles  a  =  aM  =  0  deg  and  a  -  8  deg  were  chosen,  as 
was  a  series  of  values  of  k .  Equation  (12)  was  then  used  to  calculate  a 
comparable  series  of  values  of  A,  and  then  Eq.  (48)  was  used  to  calculate 
cn-  The  circled  points  in  Figure  72  were  taken  from  the  original  tabula¬ 
tion  of  Reference  10,  which  has  already  been  discussed  in  Appendix  II.  In 
this  case,  the  agreement  between  the  tabulated  results  and  the  theory  is 
quite  good. 
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An  analytical  study  was  conducted  to  determine  if  available  unsteady  aerodynamic 
test  data  could  be  used  with  existing  helicopter  rotor  aeroelastic  and  variable  inflow 
analyses  to  predict  the  stall  flutter  response  of  a  helicopter  rotor  blade.  Incom¬ 
pressible  unsteady  aerodynamic  data  for  an  NACA  0012  airfoil  executing  pure 
sinusoidal  pitching  motions  were  employed.  To  apply  such  data  under  rotor  blade 
operating  conditions  where  multiharmonic  motions  and  velocity  variations  exist, 
the  data  were  expressed  as  functions  of  instantaneous  section  angle  of  attack  and 
its  derivatives.  In  addition,  scaling  procedures  were  developed  in  an  attempt  to 
account  for  the  effects  of  compressibility. 

Limited  application  of  the  resulting  analysis  to  define  the  aeroelastic  characteristics 
of  several  blade  designs  showed  that  significant  self-excited  torsional  oscillations  of 
the  stall  flutter  type  could  be  predicted  for  certain  combinations  of  flight  conditions 
and  blade  designs.  Control  load  correlation  studies  were  performed  with  CH-53A 
maneuvering  flight  test  data  and  with  level  flight  test  data  from  the  NH-3A  (S-61FL 
Although  the  analyses  produced  self-excited  stall  flutter  oscillations,  the  degree  of 
correlation  indicates  the  need  for  further  development.  Recommendations  for  such 
development  are  described  in  the  report. 


DMH  a  a  m  MAIACU  DO  roms  ini,  i 
OMH.ITI  PM  tIMV  Utl. 


JAM  34.  WHICH  I* 


Unclassified 


«< 


Unclassified 

Ow«lflniWi» 


■  tv  VOIKOS 


Helicopter  Rotor  Blades 
Rotor  Aeroelasticity 
Rotor  Blade  Control  Loads 
Unsteady  Aerodynamics 
Stall  Flutter 
Dynamic  Stall 


